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FOREWARD 


'After  decades  of  theoretical  development,  fog  models  have  evolved 
to  the  point  of  qualitative  agreement  with  observations.  Current 
models  are  capable  of  simulating  observed  quasi-periodic  oscillations 
in  fog  variables  as  well  as  the  various  stages  of  fog  growth.  It  is 
now  anticipated  that  continued  efforts  towards  greater  model  sophisti¬ 
cation  will  lead  to  at  least  limited  quantitative  agreement  with 
observations  and  improved  predictive  capability. 

Application  of  very  high  resolution  Landsat  digital  data  to  fog 
episodes  has  revealed  the  stratocumul us  nature  of  fog.  That  other 
satellite  imagery  has  shown  little  hint  of  this  cellular  structure  is 
due  to  the  fact  that  the  fog  elements  are  generally  smaller  than  the 
spatial  resolution  available  on  these  satellite  platforms.  The  ratio 
of  wavelength  between  cells  and  the  height  of  the  boundary  layer  is  in 
agreement  with  values  obtained  for  Benard  cells  and  longitudinal  rolls 
in  cloud  systems.  The  processes  occurring  in  fogs  are  similar,  but  on 
a  smaller  scale,  to  those  responsible  for  these  boundary  layer  cloud 
systems.  Therefore,  further  studies  of  fog  processes  are  expected  to 
provide  quantitative  information  concerning  the  processes  responsible 
for  Benard  cells,  longitudinal  rolls  and  cloud  top  entrainment 
instability. 
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1.  STATEMENT  OF  PROBLEM 


Fog  is  a  mesoscale  boundary  layer  phenomena  which  generally  forms 
under  high  pressure  and  stable  conditions  after  nocturnal  cooling 
lowers  the  near-surface  temperature  to  the  dew  point.  However,  the 
process  of  fog  formation  is  in  dispute.  Observations  in  England  by 
Roach  et  al .  (1976)  showed  that  fog  forms  only  after  wind  speeds 
decrease  to  low  values,  so  that  turbulence  is  nearly  non  existent.  On 
the  other  hand  Rodhe  (1962)  suggested  that  fog  formation  is  the  reuslt 
of  the  turbulent  mixing  of  nearly  saturated  eddies  at  different  tem¬ 
peratures.  Lala  et  al .  (1982)  found  that  freshening  of  calm  winds, 
causing  reduced  stability  and  increased  turbulence,  contributes  to  fog 
formation. 

Roach  (1976)  was  the  first  to  observe  that  fogs  often  have  quasi - 
periodic  oscillations  of  longwave  radiation,  surface  temperature, 
visibility  and  wind  speed.  Gerber  (1981)  found  that  dense  fog  was 
composed  of  highly  structured  and  sharply  defined  local  regions. 

Other  observations  by  Jiusto  and  Lala  (1983)  have  identified  five 
distinct  stages  of  fog  develoment  and  dissipation. 

Fog  prediction  models  also  have  grown  in  sophistication  but,  to 
date,  no  model  has  reproduced  observations  with  any  degree  of  satis¬ 
faction.  Mutual  interactions  of  turbulence,  radiation  and  thermody¬ 
namics  make  fog  simulations  a  very  complicated  problem.  The  choice  of 
turbulence  formulation  appears  to  be  the  single  most  important  factor 
determining  fog  development. 

The  purpose  of  the  present  study  has  been  to  develop  improved 
models  capable  of  simulating  the  qualitative  behavior  of  radiation 
fogs.  The  observations  of  Jiusto  and  Lala  at  Albany,  New  York,  have 
been  used  for  initiation  and  verification.  The  investigation  was 
directed  towards  comparison  of  predictions  made  with  a  whole  range  of 
increasingly  more  complicated  turbulence  modes.  The  goal  was  to  iden¬ 
tify  the  least  complex  model  capable  of  producing  reasonably  accurate 
predictions. 


2.  SUMMARY  OF  RESULTS 

Due  to  their  inherent  simplicity,  and  also  due  to  the  fact  that 
such  approaches  are  computationally  less  expensive  than  higher-order 
closure  techniques,  similarity-theory  turbulence  parameter! zations 
have  been  widely  used  in  the  meteorology  community.  This  was  the 
approach  taken  in  the  present  studies.  It  was  surprising  to  see  how 
well  stability-corrected  versions  of  similarity  theory  were  able  to 
produce  fog  profiles  which  were  in  qualitative  agreement  with  the 
Jiusto  and  Lala  observations. 

Six  different  turbulent  exchange  formulations  were  applied,  based 
either  upon  measurements  or  upon  si mpl i f i cat i ons  of  second-order  clo- 
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sure  approximations.  Three  different  Monin-Obukhov  stability  func¬ 
tions  were  specified,  along  with  different  parameterizations  of 
Richardson  number  and  mixing  length.  The  subset  of  results  presented 
in  Appendix  1  shows  the  range  of  variability  encountered  by  using 
these  different  formulations  for  turbulent  mixing. 

While  there  are  some  large  differences  between  the  results,  the 
similarities  are  more  important.  Note  that  the  various  models  simu¬ 
late  the  five  stages  of  fog  development  as  observed  by  Jiusto  and 
Lala.  Second,  the  quasi-periodic  behavior  of  the  fog  variables  also 
is  simulated.  Considering  the  crudeness  of  the  exchange  parameteriza¬ 
tions,  the  qualitative  agreement  with  observations  indeed  is  surprising. 
The  results  show  that  fog  formation  and  development  is  directly  corre¬ 
lated  with  the  magnitude  of  the  eddy  mixing  coefficients.  Larger  tur¬ 
bulence  generation  leads  to  more  rapid  fog  development  and  to  lower 
visibilities.  During  the  mature  fog  stage  a  series  of  fog  dissipation 
and  redevelopment  episodes  occur.  In  the  upper  portion  of  the  fog, 
radiative  cooling  causes  development  of  large  liquid  water  contents 
during  turbulently  quiet  periods.  Subsequent  destabilization  mixes 
the  upper  level  liquid  water  to  the  surface,  creating  the  observed 
periods  of  surface  fog  intensification.  A  thorough  report  of  this 
behavior  is  found  in  Appendix  1. 

The  periodicity  of  15-20  minutes  observed  in  the  various  fog 
parameters  is  not  reproduced  in  the  simulations.  Rather,  depending 
upon  the  turbulence  formulation,  periodicities  on  the  order  of 
40-60  minutes  are  predicted.  However,  the  IANDSAT  observations  offer 
a  possible  explanation.  Considering  the  typical  size  of  fog  cells, 
the  fact  that  their  typical  horizontal  aspect  ratio  is  a  value  of  two, 
and  that  typical  wind  speed  in  fogs  is  on  the  order  of  2m/s,  advection 
of  fog  cells  across  a  fixed  site  would  cause  a  variability  in  fog 
parameters  of  16-18  min.  Details  are  given  in  Appendix  2. 

A  number  of  investigators  have  studied  the  relationship  of  cloud 
element  size  and  the  height  of  the  boundary  layer.  The  ratio  of  wave¬ 
length  to  the  height  of  the  boundary  layer  in  cloud  systems  is  typi¬ 
cally  in  the  range  of  2-6.  For  fogs  we  found  values  of  2-3.  While 
the  results  are  tentative  due  to  the  small  number  of  cases  studied,  it 
is  suggested  that  processes  occurring  in  fogs  are  similar,  but  on  a 
smaller  scale,  to  t’.ose  processes  responsible  for  Benard  cloud  cells, 
longitudinal  rolls,  and  cloud  top  entrainment  instability.  It  is 
anticipated  that  improved  understanding  of  fog  processes  will  lead  to 
corresponding  improvement  in  our  understanding  of  cloud  processes. 
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(Manuscript  received  5  November  1985) 

ABSTRACT 

There  is  considerable  controversy  in  the  literature  concerning  fog  formation.  One  set  of  observations  suggests 
that  fog  forms  during  a  lull  in  turbulence,  while  another  set  of  observations  suggests  that  increased  turbulence 
leads  to  fog  formation. 

A  number  of  ftrst-order  closure  techniques  are  applied  to  numerical  simulations.  The  results  show  that  fog 
formation  and  development  is  directly  correlated  with  the  magnitude  of  the  eddy  mixing  coefficients.  Larger 
turbulence  generation  leads  to  more  rapid  fog  development  and  to  Larger  liquid  water  contents.  The  rate  at 
which  the  fog  top  grows  is  directly  related  to  the  rate  at  which  turbulence  lifts  the  inversion. 

During  the  mature  fog  stage,  a  series  of  fog  dissipation  and  redevelopment  episodes  occur.  Liquid  water 
develops  in  the  upper  regions  of  the  fog  during  the  turbulently  quiet  periods.  Subsequent  destabilization  of  the 
atmosphere  increases  turbulence  generation  and  mixes  the  upper-level  liquid  water  to  the  surface,  creating 
surface  fog  intensification.  Quasi-penodic  oscillations  in  fog  parameters  are  largest  in  the  upper  regions  of  the 
fog  and  become  progressively  damped  in  the  lower  regions  of  a  thick  fog. 

These  results  are  in  qualitative  agreement  with  the  observations  reported  by  Jiusto  and  Lala  and  support  the 
hypothesis  that  there  are  distinct  stages  of  fog  development. 


1.  Introduction 

During  the  last  half-century,  there  have  been  nu¬ 
merous  fog  field  programs.  Fog  prediction  models  also 
have  grown  in  sophistication,  but,  to  date,  no  model 
has  r< .  jduced  observations  with  any  degree  of  satis¬ 
faction.  Much  of  this  lack  of  apparent  progress  may  be 
traced  to  shortcomings  of  our  understanding  of  the 
nocturnal  boundary  layer  and  associated  parameter- 
izations  of  turbulence  generation.  Lack  of  a  clear  un¬ 
derstanding  of  how  local  factors,  such  as  moisture  and 
terrain,  contribute  to  fog  formation  also  hinders  prog¬ 
ress. 

Generally,  radiation  fog  forms  under  high  pressure 
and  stable  conditions  after  nocturnal  surface  cooling 
lowers  the  temperature  to  the  dew  point.  Observations 
in  England  by  Roach  et  al.  (1976)  show  that  fog  forms 
only  after  wind  speeds  decrease  to  low  values,  so  that 
turbulence  is  nearly  nonexistent.  Once  the  fog  forms, 
radiational  cooling  of  the  droplets  causes  a  develop¬ 
ment  and  thickening  of  the  fog. 

However,  the  process  of  fog  formation  remains  in 
dispute,  with  other  observations  suggesting  that  ces¬ 
sation  of  turbulence  may  not  be  required.  Emmons 
and  Montgomery  (1947)  and  Fleagle  (1953)  have  sug¬ 
gested  that  fog  formation  may  be  onset  merely  by  the 
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radiative  cooling  of  air  near  a  cold  surface.  On  the  other 
hand,  Taylor  (19 17)  and  Rodhe(l962)  suggest  that  fog 
formation  is  the  result  of  the  turbulent  mixing  of  nearly- 
saturated  eddies  at  different  temperatures.  Studies  of 
the  microstructure  of  fog  by  Gerber  (1981)  support 
this  hypothesis.  Indeed,  fog  observations  at  Albany. 
New  York,  by  Lala  et  al.  (1982)  suggest  that  freshening 
of  calm  winds,  causing  reduced  stability  and  increased 
turbulence,  contributes  to  fog  formation — in  apparent 
disagreement  with  those  observations  in  England  by 
Roach  et  al.  (1976).  Therefore,  at  this  point  in  time, 
there  is  considerable  uncertainty  as  to  the  processes 
responsible  for  fog  formation.  There  have  even  been 
some  reports  of  fogs  first  observed  several  meters  off 
the  ground,  suggesting  that  there  may  be  other  govern¬ 
ing  factors  responsible  for  fog  formation. 

Another  topic  of  considerable  interest  which,  as  yet. 
has  not  been  fully  explained  concerns  quasi-periodic 
oscillations  in  some  observed  parameters.  An  analysis 
of  the  observations  of  Roach  et  al.  (1976)  by  Roach 
(1976)  showed  significant  phase  relationships  between 
longwave  radiation,  surface  temperature,  and  wind 
speed  th3t  persisted  for  several  hours.  Similar  short¬ 
term  variations  of  fog  parameters  also  have  been  re¬ 
ported  by  Lala  et  al.  (1978).  Gerber  ( 198! ).  Lala  et  al. 
(1982)  and  Mason  ( 1982).  The  reponed  periods  of  os- 
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dilation  range  from  10-30  min.  with  most  observations 
in  the  range  of  15-20  min.  Lala  et  al.  (1982)  suggest 
that  nights  with  stronger  atmospheric  stability  have 
lighter  winds,  shorter  periods  of  osdllation,  and  greater 
likelihood  of  patchy  fog. 

It  is  well  known  that  although  fog  may  form  nearly 
simultaneously  over  large  regions,  there  is  considerable 
local  structure  in  fog.  Gerber  (1981)  found  that  dense 
fog  was  composed  of  highly  structured  and  sharply  de¬ 
fined  local  regions.  Roach  et  al.  (1976)  suggest  that 
variations  in  the  nature  and  slope  of  the  local  terrain 
are  at  least  partly  responsible  for  the  observed  local 
structure.  Roach  (1976)  attributes  the  quasi-periodic 
oscillations  to  gravity  waves  propagating  at  the  top  of 
the  boundary  layer.  He  establishes  a  rule-of-thumb  that 
variations  caused  by  advection  are  on  the  order  of  5- 
10  min.  while  the  longer  quasi-periodic  oscillations  are 
controlled  by  a  fluctuating  balance  between  radiation 
cooling  and  turbulent  diffusion. 

Landsat  very  high  spatial  resolution  digital  data  have 
been  analyzed  for  selected  fog  episodes.  In  each  case, 
the  fog  w  as  found  to  be  composed  of  cellular  dements 
with  the  fog  resembling  a  stratocumulus  cloud.  An 
analysis  of  the  sizes  of  these  cellular  elements  is  given 
by  Welch  and  Wielicki  (1986).  The  possibility  exists 
that  advection  of  these  cells  across  the  site  may  con¬ 
tribute  to  at  least  some  of  the  observed  variation. 
Choularton  et  al.  (1981)  found  pronounced  periodic 
fluctuations  in  liquid  water  content  with  periods  on 
the  order  of  70  s.  possibly  associated  with  convective 
motions  in  the  form  of  Benard  cells. 

To  date,  fog  models  have  not  been  successful  in  pre¬ 
dicting  quasi-periodic  oscillations.  One  purpose  of  the 
present  paper  is  to  demonstrate  that  even  simple  sim¬ 
ilarity  theory  approaches  for  the  parameterization  of 
turbulent  mixing  are  capable  of  simulating  such  be¬ 
havior. 

Due  to  their  inherent  simplicity,  and  also  due  to  the 
fact  that  such  approaches  are  computationally  less  ex¬ 
pensive  than  higher-order  closure  techniques,  similar¬ 
ity-theory  parameterizations  are  widely  used  in  the 
meteorology  community.  Indeed,  a  number  of  studies 
have  shown  that  similarity  approaches  at  least  approx¬ 
imately  simulate  the  boundary  layer  exchange  pro¬ 
cesses;  for  instance,  see  Panofsky  and  Dutton  (1984). 

While  similarity  approaches  may  provide  a  reason¬ 
able  simulation  of  turbulent  mixing  in  clear  air,  it  is 
questionable  whether  or  not  such  approaches  are  ad¬ 
equate  in  fogs.  As  a  fog  develops,  radiational  cooling 
shifts  from  the  ground  to  the  fog  top,  thereby  causing 
the  surface  inversion  to  lift.  The  establishment  of  su- 
peradiabatic  lapse  rates  near  the  ground  with  weak 
convective  motions  (Mason.  1982)  often  is  associated 
with  fog  development.  Later,  these  superadiabatic  lapse 
rates  may  be  found  near  the  fog  top  (Jiusto  and  Lala. 
1983a).  Clearly,  methods  that  calculate  Obukhov 
lengths  and  Richardson  numbers  based  upon  laser  av¬ 
erages  near  the  surface  are  incapable  of  simulating  such 


complicated  vertical  structure.  Brown  and  Roach 
(1976)  suggest  that  the  constant  flux  layer  may  only  be 
a  few  meters  deep  in  fogs,  with  patchy  regions  of  tur¬ 
bulence  above  this  layer. 

Nevertheless,  and  in  spite  of  these  limitations,  it  is 
perhaps  surprising  to  see  how  well  stability-corrected 
versions  of  similarity  theory'  are  able  to  produce  fog 
temperature  and  liquid  water  profiles  that  are  in  qual¬ 
itative  agreement  with  observations.  At  present,  there 
have  been  no  detailed  comparisons  between  higher- 
order  turbulence  formulations  and  fog  observations. 
Indeed,  due  to  the  sparsity  of  turbulence  measurements 
in  fogs,  there  is  great  uncertainty  as  to  the  prefer  pa¬ 
rameterization  of  various  coefficients  required  in  the 
higher-order  closure  formulations.  It  has  not  been 
demonstrated  that  more  sophisticated  approaches  will 
actually  produce  better  agreement  with  observations, 
even  though  this  is  certainly  expected. 

The  present  paper  applies  a  number  of  stability-cor¬ 
rected  turbulent  exchange  methods  to  simulate  fog  for¬ 
mation,  demonstrates  that  these  approaches  do  predict 
quasi-periodic  oscillations  in  fog  parameters,  and  an¬ 
alyzes  the  interaction  of  radiational  cooling  and  tur¬ 
bulent  mixing.  A  later  paper  (Part  II)  will  apply  local 
Richardson  number  approaches  and  higher-order  clo¬ 
sure  methods.  Comparison  among  the  various  ap¬ 
proaches  and  with  observations  can  then  determine 
the  cost  effectiveness  of  the  more  sophisticated  meth¬ 
ods.  If  it  should  be  demonstrated  that  these  other  ap¬ 
proaches  only  marginally  improve  the  agreement  with 
observations,  then  that  should  be  a  clear  sign  that  the 
importance  of  other  factors  such  as  wetness,  v  egetation, 
the  slope  of  the  ground,  local  topography,  and  pollution 
need  to  be  more  fully  addressed  (Mason.  1982). 

2.  Observations 

The  observations  at  Albany.  New  York,  identified 
five  distinct  stages  in  the  evolution  of  a  dense  fog  (Jiusto 
and  Lala.  1983a). 

a.  Sundown  stage 

The  surface  layer  was  observed  to  change  from  a 
lapse  condition  just  before  sunset  (1800  LST)  to  a 
strong  temperature  inversion  structure  (7°C/16  m)  bv 
2000  LST. 

b.  Conditioning  stage 

This  second  period  lasted  from  6  to  9  hours,  with 
surface  layer  cooling  decreasing  throughout  the  period 
to  about  1°C  h  '.  Cooling  was  observed  to  propagate 
to  higher  levels  as  heal  was  lost  by  downward  turbulent 
flux  to  the  cold  surface  and  radiative  fluxes  upward. 
Relative  humidity  increased  (up  to  a  height  of  16  m) 
to  values  in  excess  of  97'";  during  this  penod  Highlv 
oscillatory  behavior  in  visibility,  relative  humiditv.  and 
drop  concentrations  was  observed 
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c.  Mature  stage 

Dense  fog  showed  variable  oscillations  in  the  ob¬ 
served  variables,  but  they  became  highly  damped. 
Quasi  steady-state  conditions  prevailed  during  this 
stage.  Vertical  temperature  structure  tended  to  range 
between  isothermal  conditions  near  the  ground  to  even 
superadiabatic  lapse  rates. 

d.  Sunrise  stage 

The  dense  ground  fog  was  observed  to  become  more 
intense  shortly  after  sunrise.  Reformation  of  the  fog 
was  observed  to  occur  after  sunrise  in  other  cases. 

e.  Fog  dissipation  stage 

Radiation  fog  typically  “bums  off”  within  2-4  hours 
after  sunrise,  predominantly  due  to  heating  of  the 
ground  and  increased  turbulent  mixing  with  dry  air 
from  aloft. 

3.  The  fog  model 

The  model  simulations  reported  in  section  4  are 
compared  to  the  various  stages  of  fog  formation  listed 
in  section  2.  It  is  shown  that  the  numerical  model  pre¬ 
dicts  these  five  fog  stages. 

a  Mathematical  formulism 

1 )  Atmospheric  prognostic  equations 

Potential  temperature,  humidity,  liquid  water,  pol¬ 
lution  concentration,  and  wind  vertical  profiles  in  the 
boundary  layer,  respectively,  are  assumed  to  be  de¬ 
scribed  by  the  set  of  one-dimensional  equations  listed 
here  (Roach  et  al.,  1976;  Forkel  et  al„  1984): 


&\_  (PJPY  aF"  .  (■ 
dt  d:  l  *  d:J  pcp  d: 
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height  coordinate,  time 
potential  temperature 
humidity  mixing  ratio,  liquid  water 
air  density,  partial  density  (concentra¬ 
tion  of  dry  aerosol) 
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rate  of  condensation  per  unit  mass  of  air 
gravitational  settling  flux  of  liquid  water 
microturbulent  diffusion  coefficients  for 
fluxes  of  heat,  momentum,  humidity, 
liquid  water  and  aerosol 
horizontal  wind 
geostrophic  wind 
Coriolis  parameter 
vertical  unit  vector 
air  pressure,  standard  air  pressure 
Poisson's  constant 

net  flux  of  short-  and  longwave  radiation 
heat  capacity  of  air  at  constant  pressure 
Obukhov  length,  latent  heat 
heat  and  pollution  sources. 


2)  Diagnostic  equations 

The  present  formulism  does  not  allow  for  supersat¬ 
urated  conditions,  does  not  solve  the  droplet  growth 
equation,  and  does  not  use  aerosol  particles  as  con¬ 
densation  nuclei.  Rather,  the  McDonald  (1963)  satu¬ 
ration  adjustment  procedure  is  invoked  to  account  for 
latent  heat  release  and  liquid  water  formation.  Brown 
(1980)  shows  that  differences  in  the  predicted  fog  de¬ 
velopment  between  the  drop  growth  procedure  and  the 
saturation  adjustment  method  are  not  large. 

Following  Brown  and  Roach  (1976),  the  gravita¬ 
tional  flux  of  droplets  is  approximated  by 

G  =  Vrw,  (2) 

where  VT  is  the  average  droplet  terminal  velocity.  The 
linear  regression  model  of  Brown  and  Roach  (1976)  is 
used  to  parameterize  Kr(m  s_l)  as 

.  VT  =  0.0625vty  (3) 

for  liquid  water  contents  ranging  up  to  wy  =  0.3  g  m~3, 
as  typically  found  in  fogs. 

At  relative  humidities  below  1 00°o.  water  vapor 
condenses  on  aerosol  particles,  forming  a  thin  mist. 
This  process  is  included  by  assuming  a  Junge-tvpe  par¬ 
ticle  distribution,  using  Hand’s  (1976)  method  for  rel¬ 
ative  humidities  below  99.8%  to  compute  the  equilib¬ 
rium  radii.  For  relative  humidities  between  99.8cc  and 
100.0%.  the  size  distribution  of  the  haze  particles  is 
based  on  the  results  of  Neiburger  and  Chien  (I960). 

b.  Radiative  transfer 

The  generalized  two-stream  approximation  of 
Zdunkowski  et  al.  (1982)  is  used  to  compute  radiative 
fluxes  and  heating/cooling  rates.  The  parameterization 
of  absorption  and  extinction  coefficients  for  both  the 
solar  and  infrared  window  spectral  regions  is  given  by 
Forkel  et  al.  (1984).  Multiple  scattering  by  particles 
and  droplets  is  ignored  in  the  strong  infrared  absorption 
bands  of  water  vapor  and  carbon  dioxide.  However, 
an  extended  emissivitv  method  is  used  to  account  for 
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overlapping  absorption  effects  of  gases,  particles,  and 
droplets. 


c.  The  Brunt-  Vaisala  frequency 


A  number  of  parameterizations.  including  the  Rich¬ 
ardson  number  and  scale  length,  depend  upon  the  value 
of  the  Brunt-Vaisala  frequency,  N.  For  unsaturated  air, 
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For  saturated  air.  a  number  of  expressions  have  been 
suggested  (Dudis.  1972;  Fraser  et  a!..  1973;  Lalas  and 
Einaudi,  1974;  Durran  and  Klemp.  1982).  Lalas  and 
Einaudi  (1974)  give  the  form 


g  dqw 
1  +  qw  dz  ' 


(5) 


where  R  is  the  dry  air  gas  constant.  Fm  is  the  saturated 
adiabatic  lapse  rate,  and  qw,  the  total  water  mixing  ratio, 
is  the  sum  of  the  saturation  mixing  ratio,  qs ,  and  the 
liquid  water  mixing  ratio,  q /.  However,  the  above 
expression  can  be  difficult  to  evaluate  numerically. 
Durran  and  Klemp  (1982)  show  that  a  very  good  ap¬ 
proximation  to  Eq.  (5)  is 


1  +  {LeqJRT) 
l  +  (tLfqs/CpRT2) 


Id  lng  dq,\  _  dq^ 

\  dz  CpT  dz)  dz  ' 


(6) 


d.  Turbulence  formulations 

While  there  has  been  some  success  in  modeling  the 
unstable  boundary  layer,  an  effective  parameterization 
for  turbulent  exchange  remains  elusive  for  the  noctur¬ 
nal  stable  boundary  layer.  The  proven  and  tested  state- 
of-the-art  has  not  progressed  much  beyond  K-type  the¬ 
ories  (Blackadar,  1979).  As  there  have  been  no  serious 
efforts  to  measure  turbulent  variances  and  covariances 
under  saturated  conditions,  it  is  necessary  to  assume 
that  those  formulations  developed  for  unsaturated 
conditions  apply  equally  well  to  fogs.  It  is  assumed  that 
turbulence  is  continuous,  not  intermittent  (Kondo  et 
al„  1978),  that  there  are  no  terrain  effects  (Brost  and 
Wyngaard,  1978),  and  that  gravity  waves  are  not  pres¬ 
ent  (Finnigan  and  Einaudi,  1981). 

A  number  of  different  turbulent  exchange  formu¬ 
lations  are  applied,  each  based  either  upon  measure¬ 
ments  or  upon  simplifications  of  second-order  closure 
approximations.  The  first  two  methods  are  based  upon 
stability  corrections  to  expressions  appropriate  for 
neutral  conditions  (superscript  n): 


where  subscripts  m  and  h  refer  to  momentum  and  heat 
fluxes,  respectively,  the  are  the  Monin-Obukhov 
profile  functions,  and  K0  is  a  background  or  molecular 
diffusivity  term  below  which  the  eddy  diffusivitv  coef¬ 
ficients  are  not  allowed  to  fall  (Zhang  and  Anthes. 
1982);  Ko  =  2.4  X  10"5  m:  s-1.  The  first  of  these  for¬ 
mulations  is  due  to  Shir  (1973) 


(8) 


and  the  second  is  due  to  Blackadar  ( 1962) 
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Kmm  =  l2s  (9) 

where  h  is  the  height  of  the  boundary  layer.  !  is  the 
mixing  length,  and 


-m^)T 


(10) 


There  is  considerable  controversy  surrounding  the 
parameterization  of  both  h  and  /(Zilitenkevich.  1972; 
Delage,  1974;  Wyngaard.  1975;  Brost  and  Wyngaard. 
1978;  Blackadar.  1979;  Arva.  1981:  Andre  and  Mahrt. 
1982:  Nieuwstadt.  1984).  The  height  of  the  boundary 
layer  is  parameterized  (Brost  and  Wyngaard.  1978)  as 


h  »  d(u.Lo/f)'/2  (11) 


with  d  =  0.40.  in  agreement  with  Zilitenkevich  ( 1972, 
1975),  and  To  the  Obukhov  length.  However,  the  value 
of  o' depends  upon  the  wind  direction  relative  to  a  slope, 
the  angle  of  the  sloping  surface,  and  the  surface  cooling 
rate  (Brost  and  Wyngaard.  1978).  Furthermore,  under 
stable  conditions,  the  breakdown  from  turbulent  to 
laminar  flow  may  occur  at  very  low  levels  (Businger 
and  Arya,  1974),  with  alternating  layers  of  weakly  tur¬ 
bulent  and  laminar  flow- aloft  (Xing-Shengetal.,  1983). 

The  variable  mixing  length  is  often  expressed 
(Blackadar,  1979)  as 


I  =  k(z  +  2o)/[  1  +  k(z  +  --o)/A].  (12) 

Near  the  surface,  the  height  z  is  the  proper  scaling 
length,  but  when  z>  L,  Businger  and  Arya  (1974)  and 
Blackadar  (1979)  suggest  that  /  =  L  may  be  appropriate. 


z,  z  <  Lo 
L,  r  >  To. 


(13) 


An  alternate  expression  is  suggested  by  Brost  and 
Wyngaard  (1978)  in  the  form 


1  =  1  1 

1  z  lg 

where  lB  is  the  buoyancy  length  scale. 
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with  CB  =  1.69.  In  a  fog.  N  is  replaced  with  Nm  (Eq. 
6).  The  value  of  aw  is  from  Nieuwstadt  (1984), 

Ow~  1-4#»*[1  —  z/h]3/4.  (16) 

Equations  (13)  and  ( 1 4)  express  the  fact  that  in  a  stable 
boundary  layer,  the  turbulent  eddy  is  restricted  so  that 
turbulence  no  longer  “feels"  the  presence  of  the  surface 
(Nieuwstadt,  1984). 

The  Monin-Obukhov  stability  functions  [Eq.  (7)] 
also  need  to  be  specified.  However,  once  again  there  is 
considerable  controversy  as  to  which  expressions  are 
most  appropriate  (Blackadar,  1979).  We  consider  the 
following  three  forms. 

1 )  Businger  et  al.  (1971): 

unstable  conditions 

«m  =  (i  -  i5rr,/4 

**  =  (1 -9rr,/2  d7) 

stable  conditions 


<t>m  =  1  +  4.7f 


Based  upon  a  simplified  version  of  Wyngaard's 
( 1 975)  second-order  turbulence  for  the  stable  boundary 
layer,  Brost  and  Wyngaard  ( 1 978)  suggest  the  following 
form 


(22) 


with  A*  =  1.25Am. 

In  contrast,  Panofsky  and  Dutton  (1984)  suggest  the 
very  simple  form  for  stable  air, 

ku  ~ 

Kh  =  Km  =  -^  +  Ko.  (23) 

<Pm 

Based  upon  second-order  closure.  Blackadar  (1979) 
suggests  that  for  unstable  conditions. 

Km  =  /2i(l  -  21  Ri),/2j  (24) 

Kh  =  /2s(l  -  87  Ri),/2J 
and  for  stable  conditions. 


4>h  =  0.74  +  4.7f, 

where  f  =  :/Lq. 

2)  Carl  etal.  (1973): 

unstable  conditions 

*m  =  d  -  15fr,/3 

4>h  =  0.74(1  -  16fT,/2 
stable  conditions 


(18) 


K„  =  Km  = 
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Ri  Ric 
Ri  S*  Ri^. 


(25) 


For  both  the  Brost  and  Wyngaard  exchange  (22),  and 
the  Panofsky  and  Dutton  exchange  (23).  we  use  Eq. 
(24)  for  unstable  conditions. 


e.  The  Richardson  number 


<t>m  =  1  +  S( 


1 0.74  +  9.6f  +  29. 6f2  f  <  0.08 
[  1.2  +  6. If  f  >  0.08.  (19) 


3)  Dyer  and  Bradley  (1982): 
unstable  conditions 

<t>m  =  (1  -  28  f)' 1/4 
<t>k  =  (1  -28f)*,/2 

stable  conditions 


4>»  =  4>m=  1  +  4.7f.  (20) 


We  also  consider  a  few  other  parameterizations  for 
turbulent  mixing  in  the  steady-state,  stable  boundary 
layer.  Businger  and  Arya  (1974)  suggest  the  following 
form  for  turbulent  mixing  coefficients  of  momentum. 


ku*:  expi-Vgfc/uj) 
I  +  4.7f 


+  Ko. 


(21) 


Here  we  take  Kh  =  Km  for  the  stable  case  and  use 
Blackadar's  (Eq.  9)  exchange  with  Businger  et  al.’s 
( 1971 )  stability  (Eq.  17)  for  unstable  conditions. 


Following  Blackadar  ( 1 979),  the  Richardson  number 
Ri  is  calculated  rather  than  Obukhov  length  Lo.  The 
Richardson  number  is  expressed  for  unsaturated  con¬ 
ditions  [Eq.  (4)]  as 

Ri  =  A7s,  (26) 

and  for  saturated  conditions  (Eq.  6)  as 

Ri  =  Nm/s.  (27) 

The  similarity  function  may  be  obtained  from  Businger 
et  al.  (1971) 

for  unstable  conditions 

=  0.74f(  1  -  15f)l/:! 

(1  ~  9f)l/2 

and  for  stable  conditions 

„.  f(0.74  +  4.7f) 

(i  +  4.7f): 

However.  Panofsky  and  Dutton  (1984)  note  that  sim¬ 
pler  relationships  given  by  Businger  (1966)  and  Pan- 
dolfo  ( 1966)  differ  significantly  only  in  nearly  neutral 
conditions  and  are  not  necessarily  less  accurate  than 

Eqs.  (28)— (29). 


(28) 


(29) 
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For  unstable  air 


For  stable  air 


f  =  Ri. 


J  1  -  5  Ri  ' 

The  present  calculations  use  a  mixture  of  these  two 
approaches,  inverting  Eq.  (29)  for  stable  conditions  and 
using  a  simplified  version  of  Eq.  (28) 


Ri  =  0.74f 


for  unstable  conditions. 


/  Soil  equations 


-pKn~ — Jv-Jw  =  0,  (36) 

oz 

where  Jv  and  Jw  are  the  soil  fluxes  of  water  vapor  and 
liquid  water,  respectively.  The  saturation  vapor  pres¬ 
sure  of  water  in  soil  is 


te) 


The  equations  for  heat  and  moisture  transport  in 
the  soil  are  given  by 

(C,  +  lb) +  loo,  |  o  £  (-Ar^  +  T.jfJ-) 

-  +  s]  +  (vP  ~  n)  Yt  <33> 

Pw(l  +fl)f  +  ft^-'=  -J_(Jv  +  J~)  (34) 


C„  T,  specific  heat  of  soil  at  constant  pressure,  soil 

temperature 

/,  \T  heat  of  evaporation  of  water  in  soil,  thermal 
conductivity 

Pw,  v,  vP  density  of  water,  liquid  water  content  of  soil, 
porosity 

f  soil  moisture  potential 

fluxes  of  water  vapor  and  liquid  water  in 
soil 

hv,  hw  partial  specific  enthalpy  for  water  vapor  and 
liquid  water 

a,  b,  c  parameters  given  by  Forkel  et  al.  (1984). 

The  moisture  potential  \p  results  from  capillary  and 
adsorption  effects  and  is  a  function  of  both  temperature 
and  soil  liquid  water.  Equations  (33-34)  are  a  coupled 
pair  involving  heat  conduction,  phase  change  and 
transport.  A  complete  description  of  the  soil  model  is 
found  in  Forkel  et  al.  (1984)  and  Sievers  et  al.  (1983). 

g  Interface  conditions 

The  atmosphere  and  soil  system  are  coupled  at  the 
soil  surface  by  interface  conditions  of  heat  and  water: 

Fn~  pCpKhf  +  Xr^ 

OZ  02 

+  (-pA*^ -./„)(/,. -yfr)  =  0  (35) 


P!  =  Ps  expl 


where  Ps  is  the  saturation  vapor  pressure  of  a  flat  water 
surface  and  Rv  is  the  water  vapor  gas  constant.  Tem¬ 
perature.  liquid  water,  and  vapor  may  be  determined 
at  the  soil  interface  using  Eqs.  (35-37).  Note  that  sed¬ 
imentation  of  fog  droplets  is  assumed  to  contribute  to 
dew  formation,  but  not  to  the  soil  interface  conditions. 
The  upward-directed  water  vapor  flux  from  the  soil 
also  contributes  to  dew  formation  (Monteith,  1957). 

h.  Boundary  conditions 

The  upper  and  lower  boundary  conditions  are  taken 
as 

u  =  u.,  v  =  0  ] 


—  =  0.35°K/100  m 
02 


J-  = 


at  z  =  3  km. 


M,  top  _  q  at  first  grid  point  above  the  fog 
4/  with  relative  humidity  <  100%. 

u,v  =  0  at  z  -  z0,  roughness  height. 


T,  =  constant 


at  :  =  -1  min  soil. 


Here  g„_,  refers  to  the  mixing  ratio  one  grid  point  below 
the  upper  boundary;  ug  is  the  geostrophic  wind. 

/.  The  numerical  procedure 

The  numerical  evaluation  is  carried  out  using  an 
implicit  scheme.  For  higher  resolution  of  the  grid  near 
the  soil  surface,  the  basic  equations  are  subjected  to 
logarithmic  and  square  root  transformations  in  the  air 
and  within  the  soil,  respectively.  Details  are  given  by 
Zdunkowski  et  al.  (1976).  The  maximum  possible  time 
step  during  fog  episodes  is  determined  by  droplet  set¬ 
tling.  Integration  time  steps  of  I  s  are  required.  In  the 
absence  of  fog.  75  s  time  steps  are  used. 

j.  Initial  conditions 

The  numerical  simulations  reported  below  refer  to 
early  October  at  Albany.  New  York.  Prediction  begins 
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at  1800  (local  time  used  throughout)  and  continues 
until  1000  the  next  morning. 

The  present  study  attempts  to  simulate  the  obser¬ 
vations  of  Jiusto  and  Lala  ( 1983a.b).  Jiusto  and  Lala 
(1982).  and  Lala  et  al.  (1982)  in  October  of  1981  and 
1982.  While  a  complete  dataset  is  not  available  from 
these  reports,  a  representative  vertical  temperature 
structure  is  taken  from  these  studies  for  the  lower  1 6 
m.  No  information  is  available  above  this  height.  Sur¬ 
face  temperature  at  1800  is  taken  to  be  14°C,  decreas¬ 
ing  to  12°C  at  a  height  of  :  =  0. 1  m.  The  atmosphere 
is  assumed  to  be  isothermal  at  I2°C  up  to  a  height  of 
r  =  0.5  km.  decreasing  to  1 1°C  at  z  =  1  km  and  to 
-3°C  at  z  -  3  km.  Soil  temperature  is  assumed  to 
increase  to  15.5°C  at  a  depth  of  r  =  -5  cm,  and  then 
to  gradually  decrease  to  15°C  at  a  depth  of  r  =  1  m. 

Relative  humidity  U  (in  percent)  is  assumed  to  be 

V  =  80  up  to  a  height  of  *  =  0.5  km,  then  decreasing 
linearly  with  height  to  L  =  60  at  r  =  3  km.  Above  that 
height,  the  average  observed  autumnal  values  are  used 
for  calculating  the  radiative  transfer.  Soil  moisture  is 
assumed  to  be  ?j  =  0.16  at  the  surface,  increasing  to 

V  =  0.2  at  a  depth  of  :  =  - 1  m.  Data  for  soil  moisture 
potential  and  hydraulic  conductivity  are  taken  from 
Hillel  and  van  Bavel  ( 1976),  and  thermal  conductivity 
is  given  by  Baver  et  al.  (1972).  Geostrophic  wind  is 
taken  for  calm  conditions  as  2  m  s_l  and  the  roughness 
height  r0  is  taken  as  1  cm.  Weak  pollution  is  assumed, 
with  ground  concentration  of  5000  aerosol  particles/ 
cm3,  corresponding  to  a  visibility  of  60  km  at  U  =  50. 
The  initial  particle  distribution  is  assumed  to  decrease 
exponentially  with  height  as  exp(-r///0)  where  HD  is 
taken  as  1  km. 

4.  Results 

a.  Color  processed  display  of  results 

Figure  1  shows  a  color  processed  display  of  temper¬ 
ature.  liquid  water  content,  heating/cooling  rate,  and 
eddy  mixing.  The  horizontal  coordinate  is  local  time 
and  the  vertical  coordinate  is  height.  Even  though  the 
fog  extended  to  heights  of  200-300  m  in  the  early 
morning  hours,  only  the  lower  150  m  is  shown  in  this 
figure.  This  choice  represents  a  compromise  in  limiting 
the  coverage  of  the  upper  region  of  the  fog  in  order  to 
improve  the  resolution  near  the  ground.  Soil  temper¬ 
ature  in  the  first  meter  below  the  surface  is  included 
in  the  first  frame.  A  color  code  is  found  below  each  of 
the  individual  frames  to  convert  the  colors  to  the  ap¬ 
propriate  numerical  values.  For  instance,  yellow  rep¬ 
resents  liquid  water  content  in  the  range  0.15-0.20  g 
m  5  (second  frame).  However,  use  of  only  eight  tem¬ 
perature  classifications  ti  e.,  colors)  is  found  to  be  in¬ 
sufficient  to  fully  visualize  the  temperature  variations. 
Therefore,  a  more  detailed  color  scheme  is  used,  as 
shown  in  Frame  I. 

It  should  be  noted  that  color  image  display  contains 
a  higher  information  content  than  standard  methods 
of  presenting  results.  Therefore,  taking  time  to  become 


familiar  with  the  rainbow  color  scheme  is  advised  be¬ 
fore  proceeding  to  discussion  of  the  results.  Note  that 
the  black  and  blue  colors  generally  represent  the  small¬ 
est  values  and  the  red  and  white  colors  represent  the 
largest  values. 

b.  Businger-Arya  exchange 

This  section  discusses  the  development  of  fog  as  pre¬ 
dicted  using  the  Businger-Arya  (1974)  eddy  exchange 
formulation  [Eq.  (25)]  using  the  Businger  et  al.  ( 1971 ) 
stability  coefficients  [Eq.  (21)].  Scaling  length  is  given 
by  Eq.  (15).  and  the  height  of  the  boundary  layer  is 
calculated  from  Eq.  (11).  Comparison  is  made  with  the 
observations  of  Jiusto  and  Lala  (section  2)  to  determine 
how  well  the  model  simulates  the  observed  stages  of 
fog  development. 

1 )  Sundown  stage 

Surface  temperature  decreases  rapidly  during  the  first 
two  hours  of  simulation,  dropping  from  14°C  at  1800 
to  9.7 °C  at  1900  and  to  8°C  at  2000.  Due  to  the  large 
changes  in  temperature  in  the  lower  meter  of  the  at¬ 
mosphere  and  the  pixel  resolution  limitations  of  the 
color  imaging  system,  it  is  difficult  to  obtain  accurate 
surface  temperatures  in  Fig.  la.  Above  the  surface,  at 
heights  of  1- 10  m.  the  decrease  in  temperature  is  much 
more  gradual,  with  a  cooling  of0.5°C  h~'.  These  values 
are  much  smaller  than  the  7 °C  h~'  decrease  in  surface 
temperature  and  the  4-6°C  h_l  change  in  the  0.5-16 
m  layer  observed  by  Jiusto  and  Lala  (1983a). 

The  height  of  the  boundary  layer  (Fig.  2a)  decreases 
from  its  original  h  =  85  to  60  m  at  1900  and  then 
rapidly  drops  to  just  a  few  meters  at  2000.  During  this 
initial  stage,  radiative  cooling  of  the  surface  (Fig.  2b) 
is  dominant,  with  turbulent  transport  suppressed  due 
to  the  presence  of  the  strong  surface  inversion  (Fig.  2c). 

2)  Conditioning  stage 

During  the  next  four  hours,  surface  cooling  decreases 
by  a  relatively  uniform  1°C  h~'  to  4.2°C  at  2400  (Fig. 
la).  This  rate  of  decrease  is  similar  to  that  observed  by 
Jiusto  and  Lala.  During  this  period,  the  surface  inver¬ 
sion  increases  in  intensity  and  in  height  (Figs,  la  and 
2a).  By  2400.  the  inversion  is  6.5°C/50  m.  After  the 
collapse  of  the  height  of  the  boundary  layer  at  2000. 
the  gradual  increase  in  surface  stress  causes  a  decrease 
in  stability  (Ri  <  Ric,  Fig.  2c).  a  gradual  increase  in 
the  height  of  the  boundary  layer  (Fig.  2a).  and  increased 
turbulent  mixing  (Fig.  Id).  Increased  mixing  then  leads 
to  increased  cooling  at  higher  levels  as  heat  is  trans¬ 
ported  to  the  ground  to  be  radiated  to  space. 

The  region  below  the  top  of  the  inversion  now  begins 
to  show  rapid  increases  in  relative  humidity.  However, 
above  this  region  of  turbulent  mixing,  relative  humidity 
changes  slowly.  This  demonstrates  that  during  this 
conditioning  stage,  the  loss  of  energy  to  the  ground  by 
turbulent  transport  is  an  order  of  magnitude  more  im- 
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FIG.  2.  (a)  Height  of  the  boundary  layer,  fi;  (b)  radiation  flux  at  the  surface;  (c)  Richardson  number.  Ri:  and 
(d)  water  column  as  a  function  of  height  and  local  time  using  the  Businger-Arya  turbulence  formulation. 


portant  than  is  the  loss  of  energy  by  radiative  cooling. 
Similar  behavior  is  observed  by  Jiusto  and  Lala. 
By  2300,  relative  humidity  exceeds  99%  to  a  height 
of  10  m. 

At  2300,  the  eddy  exchange  coefficient  decreases 
slightly  in  magnitude  and  then  intensifies  (Figs.  Id  and 
2c).  It  is  at  this  point  that  the  fog  forms.  For  the  next 
45  min,  the  ground  fog  oscillates  in  height  between  10 
and  20  m,  with  a  periodicity  of  20-30  min  (Figs,  lb 
and  2d).  Liquid  water  content  averages  about  =  0.05 
g  m~3  near  the  fog  top  to  about  0.15  g  m~3  near  the 
ground.  Radiative  cooling  rates  are  strongly  correlated 
with  liquid  water  content  in  this  shallow  ground  fog, 
ranging  from  -2°C  h"1  near  fog  top  to  almost  -5°C 
h"1  near  the  ground. 

3)  Mature  fog  stage 

During  the  next  two  hours  (until  0200),  the  turbulent 
mixing  coefficients  rapidly  increase  in  intensity  and 
height  (Fig.  Id)  as  atmospheric  stability  decreases  (Fig. 
2c)  and  the  height  of  the  boundary  layer  increases  (Fig. 
2a).  The  increase  in  fog  height  and  fog  liquid  water 
content  is  directly  correlated  with  the  increase  in  tur¬ 


bulence.  This  turbulent  region  continues  to  experience 
large  values  of  heat  transport  to  the  ground.  During 
this  period,  the  region  of  maximum  radiational  cooling 
lifts  off  the  ground  into  the  upper  regions  of  the  fog. 
Radiational  cooling  rates  on  the  order  of-3°C  h'1  to 
-4°C  h"1  are  common  (Fig.  lc). 

Near  the  ground,  liquid  water  content  slowly  in¬ 
creases  to  about  0.25-0.3  g  m'3,  being  relatively  uni¬ 
form  with  height  in  the  lower  20  m.  However,  liquid 
water  is  highly  oscillatory  in  the  upper  regions  of  the 
fog,  with  fog  top  now  reaching  to  about  50  m  (Fig.  lb). 
While  there  is  some  oscillatory  behavior  in  the  liquid 
water  content  near  the  ground,  it  is  highly  damped. 

Once  fog  forms  and  reaches  a  height  of  10-20  m  at 
2400,  surface  IR  net  flux  has  dropped  by  about  a  factor 
of  2  (Fig.  2b)  and  surface  temperature  no  longer  de¬ 
creases  (Fig.  la).  At  this  point,  the  integrated  vertical 
water  column 

w'r  =  I  w  (46) 

Jo 

is  about  1-2  g  m  J,  where  zT  is  the  height  of  the  fog 
top.  At  about  wc  =  2  g  m"2  (Fig.  2d),  the  surface  ra¬ 
diative  heat  loss  to  space  is  balanced  by  downward  heat 
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transport  in  the  atmosphere  to  the  surface  and  upward  not  observed  in  the  Albany  field  measurements.  For 
heat  transport  in  the  soil.  As  the  fog  grows  thicker,  the  the  next  four  hours,  the  numerical  simulation  shows 
surface  IR  radiative  loss  continues  to  decrease  (Fig.  a  series  of  pronounced  oscillations  in  the  fog  param- 
2b)  and  the  surface  temperature  begins  to  increase  eters.  With  each  surge  in  eddy  mixing,  the  height  of 
(Fig.  1  a).  the  boundary  layer  rapidly  expands  and  contracts  ( Fig. 

The  behavior  of  the  surface  temperature  is  in  agree-  2a).  When  the  boundary  layer  height  increases  rapidly, 
ment  with  the  observations  of  Jiusto  and  Lala  at  Al-  heat  from  the  upper  portion  of  the  fog  near  the  inver- 
bany.  Figure  3  shows  a  representative  vertical  temper-  sion  top  is  transferred  downward,  and  lower  level 
ature  profile  in  the  lower  16  m  on  the  night  of  10-1 1  moisture  is  mixed  into  the  upper  regions  of  the  fog. 

October  1982.  Note  that  the  surface  temperature  in-  The  net  effect  is  the  rapid  decrease  in  liquid  water  con- 

creases  by  2-3°C  over  a  period  of  about  2.5  hours  after  tent  in  the  lower  20  m  of  the  fog  from  nearly  wt  =  0.3 
fog  forms  at  0130.  The  atmospheric  structure  becomes  g  m~3  to  about  wt  =  0.01-0.05  g  m'3.  Above  the  height 
increasingly  unstable  as  isothermal  conditions  are  ap-  of  30  m,  fog  liquid  water  briefly  increases  to  nearly 
proached  at  about  0400.  w,  =  0.4  gm'3  (Fig.  lb).  Air  temperature  in  the  lower 

While  the  predicted  temperatures  are  warmer  than  20  m  increases  by  about  0.5-0.6°C,  but  surface  tem- 

the  observed  temperatures  shown  in  Fig.  3,  the  same  perature  increases  only  by  about  0.1  °C  (Fig.  la), 

behavior  pattern  occurs  once  the  fog  forms.  The  tern-  Therefore  a  small  surface  inversion  forms  that  tends 
perature  structure  gradually  tends  toward  isothermal  to  suppress  turbulent  mixing  (Fig.  Id:  Fig.  2c). 
conditions  in  the  lower  10-20  m  (Fig.  la),  the  Rich-  The  thick  fog  aloft  radiativelv  shields  the  surface, 
ardson  number  steadily  decreases  (Fig.  2c),  and  tur-  preventing  surface  cooling.  The  surface  warms  by  0.2- 
bulent  mixing  coefficients  increase  rapidly  (Fig.  Id).  0.5°C  h_l,  while  radiative  cooling  of  the  fog  (-1.5  to 
By  0200,  a  dramatic  change  in  the  fog  occurs  that  is  - 1  °C  h_l)  causes  a  gradual  decrease  in  air  temperature. 
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Fig.  3.  Observed  vertical  profile  of  temperature  as  a  function  of  local  time  on  IO-1I  October  1982  (From  Jiusto  and  Lala.  1983a) 
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The  net  effect  is  the  rapid  breakdown  of  the  surface 
inversion  and  the  regeneration  of  turbulent  mixing. 
Coupled  with  droplet  settling,  this  increased  turbulence 
causes  the  downward  transport  of  liquid  water  from 
aloft,  with  surface  values  once  again  reaching  wy  =  0.3- 
0.4  g  m~3. 

With  another  surge  in  turbulent  mixing,  the  cycle 
repeats.  A  series  of  strong  quasi-periodic  oscillations 
are  produced  near  the  ground.  Above  30  m.  these  os¬ 
cillations  are  less  noticeable  in  the  variation  of  liquid 
water  content.  However,  these  events  act  as  a  pumping 
mechanism  to  progressively  transport  moisture  higher 
into  the  boundary  layer.  The  fog  top  no  longer  increases 
steadily  in  height,  but  its  growth  is  now  strongly  cor¬ 
related  with  the  surges  in  turbulent  mixing.  The  period 
of  oscillation  found  with  this  model  is  30-40  min.  In 
contrast,  the  observations  show  a  period  of  light  fog 
followed  by  increased  fog  intensity  at  about  sunrise. 

The  fog  continues  to  grow  to  a  height  of  about  200 
m.  in  agreement  with  the  observations.  The  liquid  water 
column  also  shows  periodicity  with  values  ranging  from 
Hr  =  10-18  g  m'2  during  this  period.  As  the  fog  grows, 
dry  air  aloft  is  mixed  down  and  gradually  decreases 
the  maximum  liquid  water  content  found  near  the 
ground.  Droplet  settling  also  contributes  to  a  decrease 
in  fog  liquid  water.  In  contrast  to  other  fog  models,  the 
maximum  liquid  water  is  not  always  found  near  fog 
top  (Fig.  lb). 

4)  Dissipation  stage 

After  0600.  the  moisture  aloft  decouples  from  the 
surface  processes  and  a  low  level  stratus  cloud  is  formed 
from  100-225  m.  with  maximum  liquid  water  content 
of  0.12  g  m~3.  The  presence  of  the  stratus  cloud  aloft, 
radiativelv  shielding  the  top  of  the  fog  layer  below,  pre¬ 
vents  further  oscillations  in  the  fog.  The  formation  of 
a  stratus  layer  aloft  is  not  reported  in  the  Albany  ob¬ 
servations. 

With  sunrise,  the  fog  rapidly  dissipates  to  a  ground 
fog  a  few  meters  thick,  with  liquid  water  contents  on 
the  order  of  wy  =  0.01-0.02  g  m~3.  The  stratus  layer 
top  lifts  to  300  m  and  the  base  lifts  to  about  200  m. 
Cloud  liquid  water  increases  to  wy  =  0.15-0.18  g  m"3 
before  gradually  decreasing. 

Rather  than  dissipating,  the  ground  fog  grows  in 
height  from  about  5  m  to  about  10  m  as  dew  and  soil 
moisture  are  evaporated  (Fig.  lb).  Liquid  water  in¬ 
creases  to  wy  =  0.05  g  m'3  at  1000  before  beginning  to 
decrease.  The  surface  fog  finally  dissipates  after  1030, 
with  high  relative  humidity  and  haze  persisting  until 
1 100  (not  shown). 

c.  Comparison  with  observations 

Sufficient  details  are  unavailable  in  the  various  pa¬ 
pers  reporting  the  Albany  fog  observations  to  make  a 
detailed  comparison  with  the  model  results.  Neverthe¬ 


less  the  model  results  do  predict  the  same  stages  of  fog 
development. 

The  model  assumes  an  initial  temperature  profile  of 
12°C  uniform  with  height  to  500  m.  increasing  to  14°C 
at  the  ground.  This  profile  is  roughly  representative  of 
the  temperatures  reported  in  the  lower  16  m  at  1800 
on  12  October  1981  and  10  October  1982.  Fog  for¬ 
mation  occurs  at  2 1 00  on  12  October  1 98 1  and  at  0 1 00 
on  10  October  1982  and  at  2230  on  29  October  1981 
(Fig.  4).  Considering  the  uncertainty  in  the  initial  con¬ 
ditions.  the  predicted  time  of  2300  for  fog  formation 
by  the  model  appears  to  be  roughly  comparable  to  ob¬ 
servations. 

Perhaps  the  most  striking  difference  between  model 
predictions  and  observations  is  the  rate  of  temperature 
decrease  during  the  “sundown"  stage.  Observations 
show  a  very  rapid  cooling  of  the  surface  layer  and  the 
eventual  development  of  a  supercooled  fog.  In  contrast, 
the  model  predicts  a  much  slower  cooling  and  warm 
fog  conditions,  with  a  minimum  surface  temperature 
of  about  4°C.  This  difference  appears  to  be  largely  due 
to  the  rate  of  IR  cooling  in  the  surface  layer.  The  fol¬ 
lowing  factors  would  lead  to  greater  surface  cooling: 
(i)  replacement  of  the  assumed  initial  temperature  pro¬ 
file  above  16  m  with  lapse  conditions:  (ii)  decrease  the 
assumed  initial  relative  humidity  profile:  (iii)  decrease 
the  assumed  aerosol  concentration:  (iv)  decrease  the 
soil  conductivity;  and  (v)  decrease  the  assumed  initial 
soil  temperature  profile. 

Figure  4  shows  liquid  water  content  measured  during 
three  different  fog  events  at  Albany.  Each  of  these  ob¬ 
servations  was  made  at  a  height  of  about  1.5  m.  The 
oscillatory  nature  of  the  fog  can  be  seen  in  these  mea¬ 
surements  with  a  period  reported  to  be  about  20  min 
(Jiusto  and  Lala.  1982:  Lala  et  al„  1982).  In  addition, 
there  seems  to  be  a  tendency  for  some  lower  frequency 
behavior  on  the  order  of  1-2  hours.  The  29-30  October 
1981  dataset  is  particularly  interesting,  showing  epi¬ 
sodes  of  fog  development  following  by  periods  of  fog 
dissipation  near  the  ground. 

The  model  predicts  that  periods  of  surface  clearing 
often  are  associated  with  increased  liquid  water  devel¬ 
opment  near  the  fog  top  (Fig.  lb).  As  the  present  ob¬ 
servations  are  only  reported  at  shelter  height,  the  pre¬ 
dicted  fog  behavior  aloft  cannot  be  verified.  However, 
when  the  detailed  Albany  case  studies  are  completed 
and  reach  the  literature,  much  more  thorough  com¬ 
parison  with  the  model  predictions  will  be  possible. 
The  following  studies  examine  to  what  degree  the  re¬ 
sults  shown  in  Fig.  I  are  sensitive  to  the  form  of  the 
eddy  exchange  formulation. 

d.  Richardson  number 

In  the  previous  calculations,  the  Richardson  number 
is  computed  over  a  height  ofr  =  1.5  m.  with  Eqs.  (29) 
and  (32)  inverted  to  obtain  L  and  f.  This  is  similar  to 
r  =  0.75  m  used  by  Brost  and  Wyngaard  ( 19"8).  How- 
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ever,  the  value  of  z  chosen  for  the  surface  layer  is 
somewhat  arbitrary  .  In  order  to  examine  the  sensitivity 
of  the  results  to  the  choice  of  z,  several  different  values 
are  used. 

As  the  value  z  of  the  surface  layer  is  increased,  the 
time  of  fog  formation  increases.  With  the  choice  of 
z  =  5  m.  fog  formation  is  delayed  nearly  two  hours, 
until  0100.  A  larger  value  of  r  leads  to  overall  greater 
atmospheric  stability  and  smaller  eddy  exchange  coef¬ 
ficients. 

The  time  of  fog  formation  is  directly  correlated  to 
the  onset  of  turbulence  generation.  No  fog  is  generated 
when  turbulence  is  suppressed.  However,  as  the  eddy 
exchange  coefficients  exceed  values  of  about  0.01  nr 
s'1,  then  fog  is  formed.  These  results  are  in  agreement 
with  the  observations  of  Jiusto  and  Lala  (1983a)  which 
show  that  increased  turbulence  and  reduced  stability 
contribute  to  fog  formation. 

The  model  predicts  that  the  inversion  lifts  off  the 
ground  with  development  of  the  fog.  Changes  in  sta¬ 
bility  are  more  rapid  with  this  z  =  5  m  formulation, 
leading  to  a  more  rapid  intensification  of  eddy  mixing. 
This  leads  to  a  fog  that  grows  in  height  more  quickly, 
and  one  with  larger  liquid  water  content.  Values  of 
w,  =  0.35-0.40  g  m  3  are  now  predicted  for  the  lower 
20-40  m. 

The  highly  oscillatory  behavior  found  in  Fig.  1  re¬ 
mains.  but  the  period  of  oscillation  lengthens  and  the 
near-surface  liquid  water  content  tend  to  be  larger. 
From  0330  to  0600,  a  series  of  oscillations  in  liquid 
water  content  occur  near  the  surface,  ranging  from  w, 
=  0.05  g  m~3  to  wy  =  0.30  g  m  3  with  a  decrease  in 
periodicity  to  —40-60  min. 

e  Stability  parameterization 

The  previous  calculations  have  been  made  using  the 
Businger  et  al.  (1971)  Monin-Obukhov  stability  func¬ 
tions,  <t>m  and  <j>h.  However,  there  is  considerable  con¬ 
troversy  as  to  which  expressions  are  most  appropriate 
for  the  nocturnal  boundary  layer.  A  rather  thorough 
discussion  is  given  by  Blackadar  (1979).  Replacement 
of  Eq.  (17)  with  the  parameterization  suggested  by  Carl 
et  al.  (1973)  given  in  Eq.  (19)  has  only  minor  impact. 
The  only  significant  difference  is  that  the  number  of 
oscillations  between  0200  and  0600  decreases  from  7 
(Fig.  lb)  to  5;  that  is.  the  period  between  oscillations 
increases. 

Replacement  of  Eq.  (17)  with  the  parameterization 
suggested  by  Dyer  and  Bradley  (1982)  given  in  Eq.  (20) 
predicts  similar  behavior  in  the  formation  and  mature 
stages.  However,  a  sudden  decrease  in  fog  top  occurs 
at  0400.  from  nearly  150  m  down  to  50  m.  The  fog 
then  grows  back  to  a  height  of  200  m  by  0600.  The 
quasi-penodic  oscillations  near  the  ground  continue 
dunng  the  period  from  0200  to  0600.  Once  again,  the 
stratus  cloud  forms  after  0600,  but  with  somewhat 
larger  liquid  water  content. 


The  conclusion  is  that  the  choice  of  stability  for¬ 
mulation  [Eqs.  ( 1 7)— (20)]  changes  the  periodicity  of 
the  oscillations  somewhat,  but  does  not  affect  the  major 
stages  of  fog  development.  The  time  of  fog  formation, 
the  development  of  a  stratus  layer  after  0600.  and  the 
time  of  dissipation  remain  relatively  unchanged. 

/  Eddy  exchange  formulations 

Using  the  Businger  et  al.  (1971)  stability  functions 
[Eq.  (17)].  the  eddy  exchange  formulations  now  are 
varied.  Liquid  water  content,  eddy  exchange  coefficient, 
and  the  heating/cooling  rate  for  the  Blackadar  (1979) 
method  are  shown  in  Fig.  5a-c.  respectively.  Liquid 
water  content  alone  is  shown  in  Fig.  5d-f  for  the  Brost 
and  Wvngaard  (1978),  Panofsky  and  Dutton  (1984). 
and  Shir  (1973)  exchange  formulations,  respectively. 
The  height  of  the  boundary  layer,  the  surface/radiation 
flux,  the  Richardson  number,  and  the  water  column 
are  shown  for  each  of  these  four  turbulent  exchange 
methods  in  Fig.  6. 

The  processes  during  the  sundown  and  conditioning 
stages  of  fog  development  occur  at  about  the  same  rate 
as  shown  in  Fig.  1  and  Fig.  2.  In  each  case,  the  height 
of  the  boundary  collapses  dunng  the  sundown  stage 
and  then  gradually  recovers  dunng  the  conditioning 
stage.  Fog  forms  between  2300  and  2400.  and  the  initial 
period  of  fog  growth  is  directly  correlated  with  the 
strength  of  the  turbulence  generation.  Cooling  rates  on 
the  order  of  -3  to  -5°C  h'1  are  common  near  the 
ground  in  this  ground  fog. 

The  largest  liquid  water  content  is  found  in  those 
regions  of  the  fog  with  the  largest  eddy  mixing  coeffi¬ 
cients.  The  height  of  the  fog  top  rises  with  increasing 
height  of  the  inversion.  During  this  stage  of  fog  for¬ 
mation,  the  fog  top  height,  height  of  the  inversion,  and 
height  to  which  turbulence  is  found  are  all  equal.  This 
is  the  height  of  the  boundary  layer  (Fig.  6).  which  is 
directly  correlated  with  the  top  of  the  fog  seen  in  the 
liquid  water  profiles  (Fig.  5). 

The  various  fog  parameters  show  oscillatory  behav¬ 
ior  near  the  surface  during  the  initial  fog  growth  stage. 
This  quasi-periodic  variation  is  most  evident  in  the 
Panofsky  and  Dutton  results  and  least  noticeable  in 
the  Shir  results.  However,  this  oscillatory  behavior 
generally  is  strongest  near  the  fog  top,  becoming  pro¬ 
gressively  damped  in  the  interior  regions  of  the  fog. 

The  Blackadar  (1962)  turbulent  exchange  formu¬ 
lation  also  has  been  examined  (not  shown  in  the  fig¬ 
ures).  A  very  shallow  (5-10  m)  ground  fog  forms  at 
2300.  which  persists  until  0300  before  growing.  The 
reason  that  this  fog  fails  to  develop  over  this  4  h  period 
is  that  turbulence  generation  remains  small  (K  <  0.01 
nr  s  ').  This  is  another  confirmation  that  fog  devel¬ 
opment  during  its  initial  growth  stage  is  directly  cor¬ 
related  with  the  strength  of  the  eddy  mixing  processes. 

In  two  cases,  using  the  Blackadar  (1979)  and  the 
Panofsky  and  Dutton  ( 1984)  eddy  exchange  methods. 
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fog  forms  aloft  at  a  height  of  10- 1 5  m  and  then  lowers 
to  the  surface.  In  both  cases,  the  fog  formation  aloft  is 
associated  with  strong  cooling  rates  at  these  levels  in 
the  absence  of  turbulent  mixing.  Additional  case  studies 
are  needed  to  determine  the  conditions  under  which 
fog  formation  aloft,  rather  than  at  the  surface,  is  likely. 

During  the  initial  fog  development  stage,  atmo¬ 
spheric  stability  decreases.  The  surface  loss  of  energy 
by  radiation  continues  to  decrease  as  the  water  column 
continues  to  increase.  The  result  is  an  increasingly  iso¬ 
thermal  temperature  structure  near  the  ground  and  a 
steadily  decreasing  Richardson  number.  The  initial  fog 
growth  stage  is  terminated  by  an  abrupt  suppression 
of  turbulence  generation  due  to  a  sharp  increase  in  the 
Richardson  number  as  a  surface  inversion  becomes 
reestablished.  Larger  decreases  in  liquid  water  content 
result  near  the  surface.  In  some  cases,  the  liquid  water 
seems  to  be  “pumped"  aloft,  but,  in  most  cases,  the 
mixing  of  dry  air  from  aloft  depletes  liquid  water 
throughout  all  levels  of  the  fog.  There  is  a  sharp  de¬ 
crease  in  value  of  the  integrated  water  column  and  a 
corresponding  increase  in  the  surface  radiation  flux. 

Each  of  these  case  studies  predicts  a  senes  of  fog 
dissipation  and  redevelopment  episodes  during  the 
mature  fog  stage.  These  episodes,  predicted  to  reoccur 
every  1-2  hours,  are  at  least  qualitatively  similar  to  the 
near-surface  liquid  water  measurements  shown  in  Fig. 

4.  The  models  suggest  that  liquid  water  is  generated  in 
the  upper  regions  of  the  fog  by  radiative  cooling.  Pe¬ 
riods  of  increased  turbulence  generation  due  to  desta¬ 
bilization  of  the  boundary  layer  mix  the  liquid  water 
from  aloft  to  the  surface. 

In  the  case  of  strong  and  prolonged  turbulence  gen¬ 
eration  from  0300  to  0400  as  predicted  in  the  Panofsky 
and  Dutton  formulation,  the  fog  lifts  completely  off 
the  ground  to  become  low  stratus.  Only  with  gradual 
decrease  in  strength  of  the  turbulence  generation  does 
this  low  cloud  lower  to  the  surface. 

In  each  of  these  four  case  studies,  fog  is  completely 
dissipated  after  sunrise.  In  most  cases,  a  low-level  stra¬ 
tus  forms  with  a  base  height  at  100-200  m.  However, 
evaporation  of  dew  and  soil  moisture  after  sunrise 
causes  a  reformation  of  a  shallow  ground  fog  of  height 
5-10.  Liquid  water  content  in  this  shallow  fog  is  pre¬ 
dicted  to  be  on  the  order  of  vty  =  0.05  g  m~3  in  most 
cases,  but  reaches  0.20  g  m-3  for  the  Shir  exchange 
method.  This  sunrise  stage  of  fog  thickening  and  re¬ 
formation  is  reported  by  Jiusto  and  Lala.  Final  dissi¬ 
pation  occurs  at  about  1030  in  all  cases. 

5.  Conclusions 

Fog  development  is  shown  to  be  very  sensitive  to 
the  formulation  of  turbulence  generation.  However, 
the  simulations  show  the  same  stages  of  fog  develop¬ 
ment  as  observed  by  Jiusto  and  Lala  (1983a).  These 
stages  are  defined  as  follows. 


a.  Sundown  stage 

The  surface  layer  changes  from  lapse  conditions  to 
a  strong  surface  inversion  over  a  period  of  about  two 
hours.  Turbulence  generation  is  suppressed.  However, 
the  predicted  decrease  in  predicted  surface  tempera¬ 
tures  is  about  half  the  rate  observed.  This  suggests  that 
the  assumed  initial  conditions  are  incorrect.  The  re¬ 
placement  of  the  initial  temperature  profile  by  lapse 
conditions,  the  decrease  of  the  initial  relative  humidity, 
and  the  decrease  of  the  initial  soil  temperature  gradient 
would  all  contribute  to  greater  surface  cooling. 

b.  Conditioning  end  fog  formation  stage 

Surface  cooling  decreases  at  about  1°C  h"1  and  re¬ 
mains  at  this  rate  for  several  hours.  Heat  is  transported 
downward  to  the  surface  by  weak  turbulence,  and  then 
is  radiated  away  to  space.  Relative  humidity  increases 
in  the  lower  5-20  m  to  near  saturated  conditions.  A 
very  weak  ground  fog  in  the  lower  1-5  m  occasionally 
develops  in  this  period.  However,  fog  development  is 
directly  coupled  with  the  destabilization  of  the  lower 
atmosphere  and  with  increases  in  turbulent  mixing. 
These  results  are  in  agreement  with  the  observations 
of  Lala  et  al.  (1982),  which  suggest  that  increased  tur¬ 
bulence  contributes  to  fog  formation.  They  are  not  in 
agreement  with  the  observations  of  Roach  et  al.  (1976). 
which  suggest  that  fog  forms  during  lulls  in  turbulence 
generation. 

Increased  turbulent  mixing  rapidly  lifts  the  surface 
inversion.  Fog  top  height  is  directly  correlated  to  the 
height  of  this  inversion.  Fog  intensity  (liquid  water 
content)  also  is  directly  correlated  to  the  strength  of 
the  turbulence.  The  region  of  maximum  cooling  lifts 
off  the  ground  and  becomes  located  in  the  upper  third 
of  the  fog.  Temperature,  liquid  water,  visibility,  heat 
and  moisture  fluxes,  and  cooling  rate  all  show  rapid 
variation  during  this  period,  especially  in  the  upper 
regions  of  the  fog. 

c.  Mature  stage 

The  dense  fog  stage  shows  a  series  of  fog  development 
and  dissipation  episodes.  The  period  of  these  episodes 
varies  greatly  with  turbulence  formulation.  The  Bus- 
inger  and  Arya  (1974)  exchange  formulation  predicts 
a  series  of  episodes  with  a  periodicity  of  40-60  min. 
However.  Blackadar’s  (1979)  exchange  formulation 
predicts  a  series  of  fog  episodes  at  the  surface  with  a 
periodicity  of  about  1-2  h.  This  periodicity  is  in  at 
least  qualitative  agreement  with  observations  (Fig.  4) 
This  model  predicts  that  liquid  water  intensifies  in  the 
upper  regions  of  the  fog  due  to  radiative  cooling.  It  is 
only  during  periods  of  atmospheric  destabilization  that 
increased  turbulence  mixes  this  upper  level  liquid  to 
the  surface. 
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Fig.  6.  Height  of  the  boundary  layer.  ft;  radiation  flux  at  the  surface;  Richardson  number.  Ri; 
and  water  column  as  a  function  of  height  and  local  time  using  four  turbulence  formulations. 
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d.  Dissipation  stage 

Observations  suggest  that  fog  often  intensifies  after 
sunrise.  The  model  results  suggest  that  this  is  caused 
by  increased  turbulence  generation  mixing  upper-level 
liquid  to  the  surface.  However,  the  models  often  predict 
the  formation  of  a  low-level  stratus  deck  that  prevents 
this  intensification  stage.  Evaporation  of  dew  and  soil 
moisture  supplies  sufficient  moisture  to  the  atmosphere 
either  to  sustain  the  fog  or  to  cause  reformation  of  the 
fog.  Final  "bum-off”  typicallv  occurs  between  1000 
and  1 100. 

The  various  first-order  closure  turbulence  models 
are  based  upon  formulations  that  are  based  only  upon 
simple  gradients.  Therefore,  it  is  surprising  to  find 
qualitative  agreement  with  observations  to  be  relatively 
good.  The  various  models  predict  quasi-periodic  os¬ 
cillations  and  longer-term  fog  dissipation  and  redevel¬ 
opment  episodes.  However,  the  numerical  simulation 
gives  no  hint  of  very  short-period  oscillations  70  s) 
observed  by  Choularton  et  al.  (1981),  which  may  be 
associated  with  convective  motions  in  the  form  of 
Benard  cells. 

These  results  are  encouraging  because  they  indicate 
that  at  least  qualitative  prediction  of  fog  events  may 
be  feasible  using  relatively  primitive  formulations  of 
turbulence  generation.  Nevertheless,  it  is  possible  that 
the  similarities  between  the  numerical  model  and  the 
Jiusto  and  Lala  observations  are  fortuitous.  For  in¬ 
stance,  variation  of  geostrophic  wind  from  2  to  3  m 
s'1  may  have  a  strong  influence  on  the  results.  In  order 
to  improve  and  to  verify  the  various  models,  a  number 
of  thorough  case  studies  is  required.  It  is  important  to 
have  a  detailed  set  of  initial  conditions  and  to  know 
the  vertical  distribution  of  turbulent  mixing,  liquid 
water,  heat  fluxes,  and  moisture  fluxes  during  these 
various  fog  development  stages.  Extension  of  the  mea¬ 
surements  to  heights  greater  than  the  16  meters  cur¬ 
rently  available  on  the  Albany  tower  would  be  most 
desirable.  The  model  results  show  that  much  of  the 
interesting  fog  development  occurs  in  the  20-50  m 
level. 
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ABSTRACT 

Landsai  Multispectral  Scanner  (MSS)  and  Thematic  Mapper  (TM)  digital  dan  are  used  to  remotely  sense 
fog  properties.  These  include  fog  cell  size  distribution,  cell  aspect  ratio  (the  ratio  of  the  length  of  the  major  and 
minor  axes  of  the  cells),  and  cell  orientation  angle  The  analysis  is  carried  out  for  four  fog  scenes,  three  high- 
inversion  radiation  fogs  in  central  California,  and  one  advection  fog  in  eastern  South  Dakota. 

Results  for  these  initial  fog  studies  indicate  that  I )  fogs  are  stratocumuius  in  naiuie.  being  composed  of 
individual  cellular  structures:  2)  the  reflectance  properties  vary  strongly  across  the  cells,  suggesting  considerable 
variation  in  liquid  water  content;  3)  fogs  often  arc  patchy,  often  revealing  surface  features  between  <og  cells;  4) 
the  ratio  of  wavelength  (A)  between  cells  and  the  height  of  the  boundary  layer  ( /;)  is  X/h  ~  2-3.  in  agreement 
with  values  obtained  for  Benard  cells  and  longitudinal  rolls  observed  in  cloud  systems;  5)  the  typical  horizontal 
aspect  ratio  of  fog  cells  is  about  a  factor  of  2;  and  6  i  observed  quasi-penodic  oscillations  of  measured  fog  variables 
may  be  caused  by  advection  of  the  cellular  structures  across  the  observational  site. 


1.  Introduction 

Fogs  often  form  almost  simultaneously  over  large 
areas  of  hundreds  of  square  kilometers.  Observation 
programs  using  towers  and  balloons  (Pilie  cl  al..  1975: 
Roach  ei  al..  1976;  Lala  ei  al..  1982;  Gerber,  1981) 
show  that  fogs  are  not  uniform,  but  rather  have  high 
frequency  variations  in  most  measured  variables. 

Roach  (1976)  was  the  first  to  observe  that  fogs  often 
have  quasi-periodic  oscillations  in  the  measured  prop¬ 
erties  such  as  temperature,  liquid  water,  visibility,  wind 
speed,  radiation  and  turbulence  generation.  Gerber 
(1981)  and  Lala  el  al.  ( 1982)  also  reported  similar  oc¬ 
currences  of  quasi-periodic  oscillations.  The  reported 
periods  of  oscillation  range  from  10-30  min,  with  most 
observations  in  the  range  of  15-20  min.  Lala  el  al. 
( 1 982)  suggested  that  nights  with  stronger  atmospheric 
stability  have  lighter  winds,  shorter  periods  of  oscilla- 
s  tion.  and  greater  likelihood  of  patchy  fog. 

Gerber  (1981)  found  that  dense  fog  was  composed 
of  highly  structured  and  sharply  defined  local  regions. 
>  Roach  el  al.  (1976)  suggested  that  variations  in  the 
nature  and  slope  of  local  terrain  are  at  least  partially 
responsible  for  the  observed  local  structure.  Roach 
(1976)  attributes  the  quasi-periodic  oscillations  to 
gravity  waves  propagating  at  the  top  of  the  boundary' 
layer.  He  also  suggested  that  variations  caused  by  ad¬ 
vection  are  on  the  order  of  5-10  min,  while  the  longer 
quasi-periodic  oscillations  may  be  controlled  by  a  fluc¬ 
tuating  balance  between  radiation  cooling  at  fog  lop 


and  turbulence  generation.  Welch  and  Ravichandran 
(1985)  came  to  a  similar  conclusion,  with  destabili¬ 
zation  of  the  boundary  layer  producing  brief  periods 
of  strong  turbulent  mixing.  However,  the  model  gen¬ 
erally  predicts  oscillations  on  the  order  of  40-60  min. 
considerably  longer  than  the  observations. 

The  present  study  focuses  upon  fog  spatial  structure 
observed  using  extremely  high  spatial  resolution  Land- 
sat  Multispectral  Scanner  (MSS)  and  Thematic  Mapper 
(TM)  digital  data.  It  is  shown  that  fogs,  at  the  time  of 
the  Landsat  overpasses,  are  composed  of  cellular  ele¬ 
ments.  That  GOES  and  AVHRR  satellite  imagery  has 
shown  little  hint  of  this  cellular  structure  is  due  to  the 
fact  that  the  fog  elements  are  generally  smaller  than 
the  spatial  resolution  available  on  these  satellite  plat¬ 
forms. 

2.  Landsat  satellite  data 

Landsat  satellites  provide  extremely  high  spatial- 
resolution  data  in  several  spectral  bands.  The  Multi¬ 
spectral  Scanner  (MSS)  instrument  has  a  nominal-spa¬ 
tial  resolution  of  57  meters  with  four  narrow  spectral 
bands  centered  at  0.55,  0.65, 0.75,  and  0.95  qm  wave¬ 
lengths.  The  Thematic  Mapper  (TM)  instrument  has 
six  shortwave  and  near-infrared  spectral  bands  centered 
at  0.48.  0.56.  0.66.  0.83.  1.65,  and  2.21  qm  with  a 
nominal  spatial  resolution  of  28.5  meters.  The  TM  also 
has  an  1 1.5  qm  spectral  band  with  1 14  meter  spatial 
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Fig.  I.  Landsat  MSS  Image  of  high-inversion  fog  near  Merced. 
CA.  at  1759  GMT  6  January  1981;  spatial  resolution  is  57  m. 


Fig.  3.  Landsat  MSS  image  of  advection  fog  near  Sioux  Falls.  SD. 
at  1640  GMT  7  November  1983;  spatial  resolution  is  57  m 


resolution.  The  Landsat  satellites  are  in  Sun-synchro¬ 
nous  orbits,  sampling  the  Earth  at  approximately  0930 
and  2130  local  time.  The  satellite  obtains  data  for  a 
1 80-km  wide  north-south  strip  along  the  ground-track. 
This  stnp  is  then  broken  up  into  170-km  X  180-km 
scenes.  Landsat  data  provides  the  spatial  resolution 
needed  to  accurately  define  the  extent  of  the  fog  episode 
and  spatial  variability  within  the  fog. 

Figure  1  shows  a  fog  episode  in  the  California  Inland 
Valley,  centered  near  the  City  of  Merced,  at  1759  GMT 
6  January  1981.  The  image  is  for  channel  7  (0.8- 1.1 


Fig.  2  GOES  image  of  fog  in  the  California  Inland  Valiev  at  1 746 
GMT  6  January  1981,  spatial  resolution  is  I  km. 


jim).  A  section  of  the  Sierra  Nevada  Mountains  is  vis¬ 
ible  in  the  upper  right-hand  corner.  The  57-meter  spa¬ 
tial  resolution  on  this  MSS  image  is  sufficient  to  show 
that  the  fog  has  considerable  structure.  However,  the 
gray-scale  used  to  process  this  image  has  been  chosen 
to  enhance  surface  features  at  the  expense  of  fog  fea¬ 
tures.  This  is  an  example  of  high-inversion  w  intertime 
radiation  fogs  which  form  when  polar  maritime  air 
stagnates  under  a  persistent  high  pressure  system  (Tay¬ 
lor,  1941;  Holets  and  Swanson.  1981). 

Figure  2  shows  a  corresponding  GOES  satellite  image 
taken  at  1746  GMT.  The  spatial  resolution  is  2  km. 
with  no  hint  that  the  fog  is  composed  of  small  cellular 
structures.  The  Landsat  scene  shown  in  Fig.  1  corre¬ 
sponds  to  a  180  X  170  km  section  outlined  approxi¬ 
mately  by  the  black  box  in  Fig.  2. 

A  second  Landsat  MSS  image  is  shown  in  Fig.  3. 
This  is  an  advection  fog  covering  eastern  South  Dakota 
and  Nebraska  at  1640  GMT  7  November  1983.  Sioux 
Falls.  South  Dakota,  in  the  upper  right-hand  corner  of 
the  image,  reported  fog  with  visibility  of  two  miles  at 
the  time  of  the  Landsat  overpass,  and  a  maximum  sur¬ 
face  wind  speed  of  14  knots  from  the  south  at  180°. 
The  body  of  water  is  the  Gavins  Point  Reservoir  on 
the  Missouri  River.  This  fog  event  was  very  w idespread, 
with  Dodge  City.  Kansas,  reporting  one-quarter  mile 
visibility  and  Abilene,  Texas,  reporting  two-mile  visi¬ 
bility. 

Figure  4  shows  two  Landsat  TM  scenes  of  the  San 
Francisco  Bay  region.  Figure  4a  shows  channel  4  (0.76- 
0.90  lira)  at  1814  GMT  31  December  1982.  and  Fig. 
4b  shows  channel  4  at  1815  GMT  1  February  1983. 
Both  cases  show  the  western  extremity  of  fog  events 
covering  the  California  Inland  Valley. 
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El(i.  4  1  jndsjl  rM  images  of  Ihe  San  Francisco  Bj\  region  ai  (ai  1814  GMT  l|  December  1482: 
and  Ihi  1815  GMT  I  f  ebruars  |48t  Spanal  resolution  is  28  5  m 


3.  Image  processing  of  the  fog  scenes 

For  each  of  the  fog  scenes,  digital  count  is  concerted 
to  nadir  reflectance  using  the  calibration  coefficients 
given  by  Robinove  ( 1  MS 2  >  for  MSS  channel  3  data  and 
by  Castle  el  al  ( 19X4)  for  the  TM  channel  4  data  Castle 
et  ai  estimate  a  calibration  accuracy  of  ?'■ .  Figure  ^ta¬ 
el  gives  senes  of  false  color  subsections  of  the  fog  event 
shown  in  Fig.  1.  I  he  color  scale  ranges  from  black  at 
zero  counts  to  white  at  97  counts.  One  color  is  assigned 
to  each  digital  count  value,  and  the  discrete  color  scale 


is  shown  at  the  lop  of  the  image.  The  scale  represents 
a  linear  range  in  pixel  reflectance  from  V~,  to  950 . 

Figure  5a  shows  a  2048  x  2048  pixel  subsection  of 
the  fog.  corresponding  to  116  km  on  a  side.  The  false 
color  processing  brings  out  the  graininess  of  the  fog. 
The  lower  right-hand  side  of  the  image  (red)  is  a  region 
of  maximum  reflectance,  indicating  that  the  fog  in  this 
region  is  optically  more  dense  This  may  be  caused  by 
thicker  fog  in  this  region,  or  fog  with  larger  liquid  water 
contents,  or  by  the  combination  The  fog  appears  to 
be  less  intense  in  the  upper  left-hand  side  of  the  image 
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and  around  the  edges  of  the  fog.  In  these  regions  one 
expects  the  highest  rate  of  fog  dissipation. 

Figure  5b  shows  a  30  km  subsection,  taken  to  the 
right  and  below  the  center  of  Fig.  5a.  The  individual 
fog  elements  become  visible  in  this  image.  Finally,  Fig. 
5c  shows  an  8  km  subsection  taken  from  the  upper 
left-hand  corner  of  Fig.  5b.  The  size  of  the  individual 
fog  elements  now  can  be  measured,  along  with  their 
orientation.  The  blue  regions  represent  clear  or  nearly- 
clear  regions  between  the  individual  fog  elements.  The 
patchiness  of  the  fog.  as  well  as  its  stratocumulus  na¬ 
ture.  is  evident  in  this  image. 

Figures  5(d-e)  give  a  corresponding  series  of  false 
color  subsections  of  the  advection  fog  shown  in  Fig.  3. 
The  color  scale  ranges  from  0  to  118  counts,  repre¬ 
senting  a  range  in  pixel  reflectance  from  3^  to  105^0. 
Figure  5d  shows  a  1 16  km  section  of  this  advection 
fog.  The  fog  is  composed  of  large  patchy  regions.  Figure 
5e  is  a  30  km  subsection  taken  from  the  center  of  Fig. 
5d.  showing  the  fog  elements  which  have  the  appear¬ 
ance  of  cloud  streets.  Figure  5f  shows  an  8  km  subsec¬ 
tion  taken  at  the  upper  boundary  just  left  of  center  of 
Fig.  5e.  The  individual  fog  streets  appear  to  be  about 
0.75  km  in  width.  The  structure  of  this  fog  is  similar 
to  that  observed  in  a  fair  weather  cumulus  cloud  street 
discussed  by  Wiclicki  and  Welch  (1986).  The  spacing 
between  cloud  streets  is  on  the  order  of  1 .5  km  for  the 
fog.  as  compared  to  about  4  km  for  the  cloud  street 
case. 

Figure  6  shows  false  color  images  of  two  San  Fran¬ 
cisco  Bay  area  fogs  (Fig.  4).  Figures  6a  and  6d  show  60 
km  subsections.  Figs.  6b  and  6e  show  30  km  subre¬ 
gions:  and  Figs.  6c  and  6f  show  8  km  subregions,  re¬ 
spectively.  The  fog  event  on  3 1  December  1982  show  n 
in  Figs.  6(a-c)  is  the  most  uniform  of  the  four  fogs 
analyzed.  The  color  scale  ranges  from  0  to  104  counts, 
representing  a  reflectance  range  from  -1*3  to  62t7. 
Variations  in  optical  depth  (i.e..  reflectance)  are  small, 
with  no  patchy  regions  in  which  the  ground  is  visible. 
An  apparent  rift  in  the  fog  is  visible  in  the  upper  right- 
hand  comer  of  Figs.  6a  and  6b.  The  cause  is  not  know  n. 
but  may  be  associated  with  shifts  in  wind  speed  and 
direction  due  to  surface  topography. 

Figure  6(d-f)  shows  a  much  more  variable  fog.  with 
some  highly  reflecting  regions  in  an  irregular  pattern. 
The  image  color  scale  ranges  from  0  to  153  counts, 
representing  a  reflectance  range  of  -1*7  to  79T.  The 
cell  structure  shown  in  Fig.  6f  is  the  smallest  of  the 
four  cases  studied.  Cell  width  is  on  the  order  of  0.3 
km  with  separations  between  cells  on  the  order  of 
0.5-0. 6  km. 

4.  Fog  cell  size  distributions 

For  57-m  Landsat  MSS  resolution  data,  even  a  0.25 
km  fog  element  will  be  composed  of  roughly  1  5  data 


pixels.  For  TM  data,  there  will  be  about  four  times  this 
number  of  pixels.  In  either  case,  there  are  sufficient 
data  pixels  to  distinguish  between  fog  edge  pixels  and 
fog  interior  pixels.  Once  the  edge  pixels  are  distin¬ 
guished,  the  individual  fog  structures  can  be  grouped 
into  individual  fog  cells.  Statistics  then  can  be  gener¬ 
ated,  grouped  into  fog  cell  size  classes.  A  thorough  ex¬ 
amination  of  this  method  is  given  bv  Wielicki  and 
Welch  (1986). 

The  procedure  is  as  follows.  First,  a  fog  threshold 
digital  reflectance  count.  \’c ,  is  set.  Each  pixel  with  a 
digital  count  greater  than  IV  is  flagged  as  a  fog  pixel. 
Second,  if  a  fog  pixel  has  at  least  one  of  its  sides  adja¬ 
cent  to  a  clear  pixel,  then  the  pixel  is  flagged  as  an  edge 
pixel.  Third,  fog  edge  and  fog  interior  pixels  are  grouped 
into  indiv  idual  fog  cells.  A  fog  element  is  determined 
to  be  unique  when  it  has  no  edge  pixels  adjacent  to 
another  fog  cell. 

The  preceding  analysis  is  accomplished  in  a  single 
pass  through  the  Landsat  satellite  data,  starting  at  the 
top  of  the  image  and  examining  three  scan  lines  at  a 
time  while  moving  from  top  to  bottom  through  the 
image.  An  example  of  this  process  is  given  by  Wielicki 
and  Welch  ( 1986).  Fog  cell  area  is  determined  as  the 
number  of  fog  pixels  multiplied  by  the  pixel  area  (57 
X  57  m  for  MSS  data  and  28.5  X  28.5  m  for  TM  data). 
Effective  cell  diameter  is  defined  as  a  circle  of  area  equal 
to  the  area  of  the  fog  element. 

All  fog  cells  in  the  image  are  grouped  into  a  finite 
number  of  size  classes  by  effective  fog  element  diameter 
D.  Because  the  number  of  fog  elements  tends  to  de¬ 
crease  exponentially  with  increasing  fog  cellular  size, 
intervals  are  chosen  in  equal  steps  of  0.4055  in  InD. 
For  this  step  size,  each  fog  size  boundary  increases  by 
a  factor  of  1.5  (i.e..  0.059.  0.088.  0.132.  0.198.  0.296. 
0.444.0.667.  1.00  km..  .  .).  This  selection  of  fog  cell 
size  classes  maintains  a  large  number  of  fog  elements 
in  each  size  class  for  statistical  reliability. 

Figure  7  shows  the  fog  cellular  size  distribution  for 
each  of  the  four  fog  case  studies.  In  order  to  best  de¬ 
lineate  cloud  cellular  structure  within  the  cloud  field, 
thresholds  1 were  chosen  so  that  one-half  of  the  fog 
area  was  brighter  than  the  threshold  l  ] .  The  analy  zed 
area  is  the  same  2048  x  2048  pixel  area  shown  in  Figs. 
5a.  5d.  6a.  and  6d.  The  fog  cell  size  distribution  is  given 
as  the  number  of  fog  elements  per  1  km:  surface  area 
per  1  km  fog  cell  diameter  range.  All  four  scenes  give 
similar  numbers  of  fog  elements  below  0.5  km  in  di¬ 
ameter.  The  South  Dakota  advection  fog  (Fig.  3)  and 
the  San  Francisco  Bay  area  fog  (Fig.  4b)  have  the 
smallest  number  of  fog  cells  in  the  intermediate  size 
classes  (0.5-2. 5  km),  while  the  patchy  central  valley 
fog  (Fig.  I )  has  the  largest  number  of  cells  in  this  size 
range.  The  more  uniform  fog  (Fig.  4a)  has  no  cell  sizes 
greater  than  14  km  in  effective  diameter.  In  contrast. 
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the  centra!  California  fog  has  cells  as  large  as  22  km. 
and  the  advection  fog  cells  are  as  large  as  27  km.  Nev¬ 
ertheless.  Fig.  7  shows  that  the  different  fogs  arc  sur¬ 
prisingly  similar  in  the  size  distribution  of  the  cell 
structures. 

Figure  8a  shows  the  fog  cell  horizontal  aspect  ratio 
as  a  function  of  cell  size  for  all  size  classes  containing 
four  or  more  fog  cells.  Aspect  ratio  is  defined  as  the 
ratio  of  the  length  of  the  major  and  minor  axes  for 
each  cell.  For  a  circular  element,  aspect  ratio  is  unity. 
Figure  8a  shows  that  the  smallest  cells  have  an  average 
aspect  ratio  of  two  for  all  four  case  studies.  The  aspect 
ratio  increases  with  increasing  cell  size,  reaching  a  value 
of  about  ten  for  the  largest  cells.  The  standard  deviation 
of  aspect  ratio  within  a  fog  cell  size  class  is  about  one 
for  the  smallest  clouds.  The  standard  deviation  increa¬ 
ses  to  about  six  for  cells  with  equivalent  diameters  of 
about  10  km.  The  advection  fog  (Fig.  3)  and  the  San 
Francisco  Bay  area  fog  (Fig.  4b)  show  large  mesoscale 
variability  and  have  larger  standard  deviations  of  aspect 
ratio.  Nevertheless,  there  is  a  clear  trend  towards  larger 
aspect  ratios  for  the  larger  cells  for  all  four  fog  cases. 
Note  that  aspect  ratios  shown  arc  for  fog  cells  with 
diameters  greater  than  twice  the  Landsat  pixel  reso¬ 
lution. 

Figure  8b  shows  the  orientation  angle  of  the  cells  as 
a  function  of  cell  size.  The  orientation  angle  is  defined 
as  0°  for  elements  aligned  vertically  in  the  image,  in¬ 
creasing  clockwise.  The  Landsat  image'  is  aligned  at 
about  10°  from  north;  therefore  10°  should  be  added 
to  each  value  shown  in  Fig.  8b  to  convert  to  earth  ref¬ 


erence.  The  smaller  elements  of  the  fogs  are  aligned  in 
a  NE-SW  orientation  at  about  40°  from  north.  Figures 
5c  and  6f  show  this  orientation.  On  the  other  hand, 
the  advection  fog  has  orientation  angles  (Fig.  8b)  of 
about  20°  from  north  for  cells  with  equivalent  diam¬ 
eters  of  0.5-2  km.  Cells  of  this  size  range  are  conspic¬ 
uous  in  Fig.  5f.  and  are  aligned  in  a  more  nearly  N-S 
orientation  than  are  the  smaller  cells.  Note  that  the 
largest  cells  in  each  fog  tend  to  prefer  a  more  E-W 
orientation.  The  fogs  in  Figs.  3  and  4a  gave  large  stan¬ 
dard  deviations  of  fog  cell  orientation  angle,  presum¬ 
ably  due  to  the  mesoscale  variability  seen  in  Fig.  3  and 
the  lack  of  well-defined  fog  streets  in  Fig.  4a. 

5.  Discussion 

A  number  of  investigations  during  the  last  five  de¬ 
cades  have  focused  upon  the  generation  of  Benard  cells 
and  longitudinal  rolls  (Ackerman.  1967;  kamiko  and 
Okano.  1971;  Kuettner,  197 1 ;  Agee  cl  al .  I973;krish- 
namurti.  I975;LeMone.  1976;  Agee  and  Lomax.  1978; 
Brown.  1980;  Agee.  1984;  Atkinson.  1981).  If  wind 
direction  and  wind  speed  vary  slowly  w  ith  height,  then 
Vinnichenko  ei  at.  (1973)  point  out  that  two-dimen¬ 
sional  cells  typically  are  formed  with  longitudinal  bands 
elongated  along  the  wind  direction. 

Wind  direction  at  Sioux  Falls.  South  Dakota,  on  7 
November  1983  was  reported  to  be  170°  near  the  sur¬ 
face.  Due  to  surface  friction,  one  expects  a  veering  of 
the  wind  of  50°-60°  under  stable  conditions.  Such  val¬ 
ues,  when  added  to  the  reported  value  of  170°.  give 
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FIG.  8.  (a)  Fog  cell  aspec!  ralio  vs  cell  equivalent  diameter,  and  (b)  cell  orientation  angle 
with  respect  to  the  top  of  the  image  vs  cell  equivalent  diameter. 

wind  directions  aloft  of  220°-230°.  Assuming  that  the  length  (X)  scales  of  organized  circulations,  with  X 

fog  elements  align  along  the  wind  direction,  values  of  =  0.3-0.5  km  due  to  the  turbulent  plumes  and  1.3- 

220°-230°  are  consistent  with  the  statistical  evaluation  1 .7  km  due  to  inflection  point  instability  (Brown.  1 972. 

shown  in  Fig.  8b.  Wind  speeds  have  not  been  evaluated  1980).  Theoretical  studies  suggest  the  ratio  X/h  in  the 

for  the  California  fogs.  range  of  1.5-4  (Kelley,  1983). 

A  number  of  investigators  have  studied  the  relation-  Under  stable  conditions  such  as  found  in  fog  epi- 
ship  of  cloud  element  size  and  the  height  of  the  bound-  sodes,  h  is  the  inversion  height.  Holets  and  Swanson 

ary  layer.  Applying  linear  stability  analysis,  Pellew  and  ( 1981)  report  observations  made  in  the  high-inversion 

Southwell  (1940)  found  that  circular  cells  were  char-  fog  episodes  of  the  California  central  valley.  They  report 
acterized  by  D  =  2. 16  /;,  where  D  is  the  cell  diameter  maximum  fog  layer  heights  of  about  500  m  and  min- 

and  h  is  the  height  of  the  boundary  layer.  Walter  and  imum  heights  of  about  200  m.  Therefore  a  value  of 

Overland  (1984)  observ  ed  fog  longitudinal  rolls  from  300  m  for  h  during  the  time  of  the  Landsat  ov  erpasses 

satellite  imagery .  They  delineated  a  hierarchy  of  wave-  seems  reasonable.  Typical  wavelengths  in  the  four  fog 
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cases  range  from  X  =  0.5-0. 6  km  for  the  "uniform” 
fog  case  (Fig.  6a-c)  to  X  =  1  km  for  the  other  fogs. 
Therefore  we  estimate  that  values  of  \/h  =  2-3  are 
representative  of  the  present  fog  cases.  Since  these  val¬ 
ues  are  consistent  with  values  found  by  Kelley  (1983). 
we  are  led  to  the  hypothesis  that  processes  occurring 
in  fogs  are  similar,  but  on  a  smaller  scale,  to  those 
processes  responsible  for  Benard  cloud  cells,  longitu¬ 
dinal  rolls,  topographically  induced  gravity  waves,  and 
cloud  top  entrainment  instability. 

Once  the  height  of  the  boundary  layer  is  know  n,  a 
first-order  estimate  of  the  vertical  profile  of  the  tur¬ 
bulent  exchange  coefficient  (A)  can  be  made  (Brost 
and  Wyngaard.  1978). 

A  =  kit, -( 1  ~  :/h)'  7</>  for  .v  =£  h  ( 1 ) 

where  k  is  von  Karman’s  constant  (k  =  0.4)  and  u,  is 
the  friction  velocity  .  Holets  and  Swanson  (1981)  report 
wind  speeds  typically  on  the  order  of  2  m  s"1  at  a  height 
of  10  m.  A  representative  value  for  the  friction  velocity 
is  estimated  to  be  t<*  =  0.1  ms'1  and  the  value  of  It 
is  taken  to  be  300  m.  Function  0  is  a  stability  correction 
factor  which  is  usually  parameterized  in  terms  of  the 
Richardson  number  ( Ri )  or  the  Obukhov  height. 

On  the  basis  of  fog  modeling,  a  value  of  Ri  =  0.1 
is  adopted  representing  the  fact  that  the  high  inversion 
fogs  exist  under  conditions  of  strong  inversions.  For 
such  stable  conditions 

0=1+  5Ri,  (2) 

or  0  =  1.5.  For  these  values.  Eq.  ( 1 )  suggests  maximum 
turbulent  mixing  coefficients  on  the  order  of  1  nr  s'1. 
Such  values  are  consistent  with  model  predictions 
(Forkel  el  ai,  1984;  Welch  and  Ravichandran,  1985) 
during  the  morning  hours  after  sunrise,  suggesting  rel¬ 
atively  isothermal  vertical  profiles  and  uniform  liquid 
water  vertical  profiles  at  least  in  the  lower  100  m  of 
the  fog.  At  present,  these  conclusions  cannot  be  verified, 
as  field  programs  have  not  emphasized  turbulence 
measurements  in  fogs.  These  large  eddy  diffusion  rates 
also  suggest  that  mixing  with  the  drier  environmental 
air  surrounding  the  fog  may  be  an  important  process 
for  dissipation.  Solar  radiation  heats  the  surrounding 
surfaces,  causing  a  circulation  pattern  to  be  set  up 
analogous  to  a  land-sea  breeze  (Gustafson  and  Was- 
serman,  1 976).  Time  lapse  satellite  images  confirm  that 
fogs  tend  to  dissipate  from  the  edges  inward  (Gurka. 
1978). 

Landsat  images  of  two  fog  events  are  show  n  in  Fig. 
9,  shortly  before  dissipation.  Figure  9a  shows  the  Sac¬ 
ramento  region  on  4  March  1979  which  had  reported 
0.25  mile  visibility  at  0700  local  time;  the  airport  re¬ 
ported  5  mile  visibility  in  haze  at  the  time  of  the  Land- 
sat  overpass.  Figure  9b  shows  the  Merced  region  on  1 1 
February  1981  with  reported  light  fog  at  dawn;  this  is 
the  same  region  as  shown  in  Fig.  I .  These  are  the  final 
remnants  of  fog  episodes  which  earlier  covered  the 
California  Inland  Valley.  The  cell  size  is  large,  indi¬ 


FlC.  9.  Landsat  MSS  images  of  fog  episodes  showing  dissipation  rem¬ 
nants  in  the  C  alifornia  Inland  Valley  lai  near  Sacramento  on  4  March 
1979.  and  (hi  near  Merced  on  1 1  February  1481 

eating  large  turbulent  mixing  ratios.  This  may  be  sup¬ 
plementary  evidence  supporting  the  hypothesis  that  the 
fog  dissipates  from  its  edges  inward. 

Finally,  we  consider  the  questions  of  the  quasi- 
pcriodic  oscillations  observed  in  the  measurements  of 
Roach  cl  al.  (1976).  Gerber  (1981).  and  L.ala  cl  a! 
(1982).  First,  assume  that  the  larger  fog  elements  typ¬ 
ically  are  separated  by  distances  on  the  order  of  X 
=  1  km.  as  seen  in  the  Landsat  images.  Second,  take 
the  wind  speed  I u)  as  2  m  s  '.  a  value  typical  of  high 
pressure  under  which  log  forms  I  he  resulting  time,  i 
=  X/f.  is  8-9  mtn.  err  about  half  the  value  reported  for 
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the  quasi-periodic  oscillations  (e,g..  Roach.  1976). 
However,  note  that  the  wavelength  X  is  measured  per¬ 
pendicular  to  the  wind  speed,  and  that  fog  elements 
typically  align  along  the  wind  direction.  Since  the  typ¬ 
ical  aspect  ratio  for  fog  cellular  elements  is  a  value  of 
2,  one  would  expect  that  the  typical  distance  between 
cellular  elements  in  the  wind  direction  is  about  2X. 
The  corresponding  time  is  then  t  =  2X/r,  or  16-18 
min.  This  value  is  in  close  agreement  with  the  observed 
quasi-periodic  oscillations  of  measured  fog  parameters. 
This  suggests  that  the  observed  quasi-penodic  oscilla¬ 
tions  may  be  a  result  of  fog  cells  advected  across  the 
observational  site. 

6.  Conclusions 

The  present  results  show  that  high  resolution  Landsat 
data  is  a  valuable  tool  in  the  investigation  of  fog  prop¬ 
erties.  The  Landsat  data  prov  ide  information  concern¬ 
ing  regional  fog  coverage,  the  patchiness  of  fogs,  cell 
size  distribution,  and  cell  aspect  ratio  and  orientation 
angle.  The  major  conclusions  of  this  preliminary  study 
are 

1 )  For  the  four  case  studies.  Landsat  imagery  shows 
fogs  to  be  composed  of  individual  cellular  elements 
and  to  have  an  appearance  similar  to  stratocumulus 
fields. 

2)  Most  cellular  fog  elements  show  considerable 
variation  in  reflectance  properties  (i.e..  optical  depth) 
due  to  variations  in  liquid  water  content  and/or  fog 
thickness.  However,  infrared  TM  images  do  not  show 
much  grainy  texture,  indicating  that  fog  thickness  is 
relative!)  constant.  Therefore  the  variations  of  reflec¬ 
tance  values  must  be  primarily  associated  with  varia¬ 
tions  in  liquid  water  content.  The  patchy  nature  of  fog 
is  clearly  evident  in  the  image  processed  scenes. 

3)  The  number  of  cells  decreases  exponentially  with 
increasing  cell  diameter,  a  characteristic  also  observed 
in  cumulus  cloud  studies  (Plank,  1969;  Wielicki  and 
Welch.  1986).  The  size  distribution  is  remarkably  sim¬ 
ilar  for  both  radiation  fogs  and  advection  fogs,  with  no 
clear  differentiation  between  the  two  types.  However, 
the  advection  fog  shows  a  stronger  tendency  for  the 
fog  elements  to  form  cloud  streets. 

4)  The  wavelength  (X)  between  the  larger  fog  ele¬ 

ments  is  on  the  order  of  I  km  for  three  of  the  case 
studies,  but  with  values  as  low  as  0.5  km  for  the  more 
"uniform”  fog  case.  Using  a  value  of  300  m  as  an  es¬ 
timate  for  the  height  of  the  boundary  layer  in  high 
inversion  fogs  (Holets  and  Swanson,  1981 ).  we  obtain 
the  ratio  \/h  2-3,  in  agreement  with  studies  of  cloud 

systems  (Kelley.  1983).  These  results  suggest  that  the 
processes  responsible  for  Benard  cells  and  longitudinal 
rolls  observed  in  clouds  probably  are  also  responsible 
for  such  structure  observed  in  the  Landsat  fog  images. 
Close  to  the  surface,  friction  limits  the  scale  of  turbulent 
mixing,  thereby  producing  cellular  elements  of  smaller 


size.  Apply  ing  arguments  for  gravity  waves  to  cellular 
convection  (Vinnichenko  ei  al,  1973),  it  appears  that 
the  shorter  the  wave  (the  distance  between  cellular  ele¬ 
ments),  the  less  stable  is  that  cellular  element. 

5)  The  results  suggest  for  conditions  typical  of  ra¬ 
diation  fogs,  that  the  quasi-periodic  oscillations  mea¬ 
sured  in  fogs  may  be  the  result  of  fog  cells  advected 
across  the  observational  site.  The  period  of  16-18  mm 
is  in  close  agreement  with  the  period  of  15-20  min 
typically  found  in  fogs.  However,  the  present  prelim¬ 
inary'  study  includes  only  four  fog  case  studies  for  only 
two  different  regions.  Considerably  more  effort  needs 
to  be  expended  to  determine  how  representative  these 
results  are  for  fogs  in  different  regions.  Landsat  obser¬ 
vations  also  are  limited  to  overpasses  in  the  morning 
hours.  Nocturnal  observations  should  be  made  to  de¬ 
termine  if  the  present  results  are  representative  of  con¬ 
ditions  before  sunrise. 
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